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ABSTRACT 

We study the spectral energy distributions and evolution of a large sample of optically selected 
quasars from the Sloan Digital Sky Survey (SDSS) that were observed in 323 Chandra images analyzed 
by the Chandra Multiwavelength Project. Our highest-confidence matched sample includes 1135 X- 
ray detected quasars in the redshift range 0.2 < z < 5.4, representing some 36 Msec of effective 
exposure. We provide catalogs of QSO properties, and describe our novel method of calculating 
X-ray flux upper limits and effective sky coverage. Spectroscopic redshifts are available for about 
1/3 of the detected sample; elsewhere, redshifts are estimated photometrically. We detect 56 QSOs 
with redshift z > 3, substantially expanding the known sample. We find no evidence for evolution 
out to z ~5 for either the X-ray photon index P or for the ratio of optical/UV to X-ray flux aox- 
About 10% of detected QSOs show best- fit intrinsic absorbing columns > 10^ cm~^, but the fraction 
might reach '^1/3 if most non-detections are absorbed. We confirm a significant correlation between 
aox and optical luminosity, but it flattens or disappears for fainter {Mb^ — 23) AGN alone. We 
report significant hardening of P both towards higher X-ray luminosity, and for relatively X-ray loud 
quasars. These trends may represent a relative increase in non-thermal X-ray emission, and our 
findings thereby strengthen analogies between Galactic black hole binaries and AGN. For uniformly- 
selected subsamples of narrow line Seyfert Is and narrow absorption line QSOs, we find no evidence 
for unusual distributions of either aox or P. 

Subject headings: galaxies: active ~ quasars: general - quasars: absorption lines - surveys - X-rays 



I. INTRODUCTION 

Interest in the properties of active galaxies and their 
evolution has recently intensified because of deep connec- 
tions being revealed between supermassive black holes 
(SMBH) and galaxy evolution, such as the relationship 
between the mass of galaxy spher oids and the SMBH 
they host (the Mbh — o' con nection; iFerrarese fc MerrifrB 
120001: iGebhardt et al.|[2000l ). A feedback paradigm could 
account for this correlation, whereby winds from active 
galactic nuclei (AGN) moderate the SM BH growth by 
truncating that of their host galaxies fe.g.. lGranato etaLl 
[200l . Feedback models may explain the correspon- 
dence between the local mass density of SMBH and the 
luminosity density produced by high redsh ift quasars 
(|Yu fc Tremaind 120021: iHopkins et aUl2006bD as well as 
the 'cosmic downsizing' (decrease in the space density 
of luminous AGN) seen in AGN luminosity functions 
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jBarger et al.l 120051 : iHasingerl 120051 : IScannapieco et all 
200^ ) If quasar activity is indu ced by massive mergers 
fe.g.. lWvithe fc Loebll200a [2005^ . then theji gsaw puzzle 
now assembling may merge smoothly with cosmological 
models of hierarchical structure formation. 

Many, if not most, of the accreting SMBH in the 
Universe may be obscured by gas and dust in the 
circumnuclear region, or in the extended host galaxy. 
The obscure d fraction rnay depend on both luminosity 
and redshift (lUeda^aIll2'003l:lBrandt fc HasiMall2005i : 
iLa Franca et al.|[2005f ). and is indeed likely to evolv e on 
grounds both t heoretical (e.g.. [Hopkins et al.ll2006al) and 
obser vational (iTreister fc Urrvl 120061 : iBallantvne et al.l 
l2006f ). Such evolution seems to be required for AGN pop- 
ulations to compose the obs erved spectrum o f the cosmic 
X-ray background (CXRB; iGilh et al.|[2007[) . However, 
a full census of all SMBH remains observationally chal- 
lenging, since some are heavily obscured, or accreting 
at very low rates, below the sensitivity limits of current 
telescopes even at low-redshifts. 

AGN unification models explain many of the observed 
differences in the spectral energy distributions (SEDs) 
of AGN as being due to the line-of-sight effects of 
anisotro pic distributions of obscuring material near the 
SMBH () Antonuccil 1 1 993f) . The intrinsic number ratio of 
obscured to unobscured AGN may evolve, and is almost 
certainly a function of luminosity. Indeed, the ratio in the 
Seyfert (low-) luminosity regime is currently estimated to 
be ~4, whereas for the Q SO (high-) lumin osity regime, 
it may be closer to unity (jCilli et al ll2007D . 

Astronomers, like most people, usually look where 
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they can see. Type 1 quasars are the easiest AGN 
to find in large numbers via either spectroscopic or 
color selection because of their broad emission lines 
(FWHM>,1000 km s~^ ) and generally blue continuum 
slopes and brighter magnitudes. Most of these show 
few other signs of obscuration such as infrared excess or 
weakened X-ray emission. Large samples of Type 1 QSOs 
- as the brightest high redshift objects - have served to 
probe intervening galaxies, clusters, and the intergalactic 
material (IGM) along the line of sight, right to the epoch 
of reionization. Because observed SEDs are thought to 
be less affected by obscuration and therefore more repre- 
sentative of intrinsic accretion physics, the evolution of 
this Type 1 sample is of interest as well. 

The SEDs and clustering properties of Type 1 QSOs 
have been studied in increasing detail, probing farther 
into the Universe and wider across the sky and the 
electromagnetic spectrum. SEDs including mid-infrared 
photometry from Spitzer were compiled and character- 
ized recently for 259 qua sars from the Sloan Digital 
Sky Survey (SDSS) by Richards et all i2006H ). These 
studies are most useful for calculating bolometric lumi- 
nosities and ii'-corrections towards understanding the 
energetics of the accretion process, and its evolution 
across cosmic time. Studies that target SDSS QSO pairs 
with small separations find a significant clustering ex- 
cess on small scales ( <40 kpc/i) of varying strengths (e.g., 
iHennawi et all l2006t iMvers et all l2007l l2008f) . which 
could be due to mutual triggering, or might simply re- 
sult from the locally overdense environments in which 
quasars form (Hopkins et al.. . 2008.1. 

Optical luminosity f unction studies of o p tically 
selected quasars (e.g., [Richards et all l2006al 120051: 



Groom et al.ll2004) date back decades Ce.g.. iBovle et alj 
1988() . Accurate luminosity functions are needed to 



trace the accretion history of SMBH and to contrast the 
buildup of SMBH with the growth of galaxy spheroids. 
An increasing number of optical surveys not only select 
AGN photometrically, but also determine fairly reliable 
photometric redshifts for them. These samples stand to 
vastly improve the available statistical reliability and the 
resolution available in the luminosity/redshift plane. 

X-ray observations have been found to efficiently se- 
lect AGN of ma ny varieties, and at higher surface den- 
sities than ever ( Hasingerl l2005f) . Independent of the 
classical AGN optical emission line criteria. X-rays are 
a primary signature of accretion onto a massive com- 
pact object, and the observed X-ray bandpass corre- 
sponds at higher redshifts to restframe energies capa- 
ble of penetrating larger intrinsic columns of gas and 
dust.* Even for Type 1 QSOs, X-ray observations have 
revealed new c onnections (e.g., between black holes from 
10 to 10^ Mg: lMaccarone et al.l[200llFalcke et al.ll2004( l 
and new physical insights such as the p ossibility of ubiq- 
uitou s powerful relativistic outflows (|Middleton et al.l 
|2007|) or of relativistically broadened fluor escent FeKa 
emiss i on (see, e.g., the review articles by iFabian et al.l 
120001 [Reynolds fc NowakI l2003h . Quasars with broad 
absorption fines (BALQSOs) blueward of their UV 
emiss ion lines turn ou t to be high l y absorbed in X- 
rays (iGreen et all [19951 : IC^ree^ fTOM loTem et al.l[200ll : 



° The observed-frame, effective absorbing column is A^jj' 
Nu/il + z)'^-'^ l|Wilman fc Fabian |[T999h . 



[Gallagher et al.[ [2002( ). Quasars that are radio loud 
are 2-3 tim es brighter in X-rays for the same opti- 
cal mag nitude (Za morani et al.|[198lMWorrall et al.[[1987l 
[Shen et al. 2006), and also m ay have har der X-ray spec- 
tra (e .g., IShastri et al.[[199"3l . although see [Galbiati et al.l 

\mm . 

The high sensitivity and spatial resolution of Chandra 
and XMM-Newton open other avenues for exploration 
of quasars and their environments. Clusters of galaxies 
have been discov ered in t he vicinity of, or along the sight- 
lines to quasars (jGreen et al. 20051: [Siemigin owska et al.l 
l2005h . Lensed quasars have now been spatially resolved 
in X-rays, unexpectedly showing sig nificantly different 
flux ratios than at other wavebands ([Green et al.l [2002[ : 
[Blackburne et al.|[2006[ : [Lamer et al.| [2006'). 

The expected evolution in the environment, accretion 
rates, and masses of SMBH in AGN should correspond 
to observable evolution in their SEDs. The two most 
common X-ray measurements used are the X-ray power- 
law photon index and the X-ray-to-optical spectral 
slope, aoa;-^*^- Many of the apparent correlations have 
been challenged as being artifacts of selection or the by- 
products of small, heterogeneous samples, which impedes 
progress in our understanding of quasar physics and evo- 
lution. 

The archives of current X-ray imaging observatories 
like Chandra and ^MM-Newton are growing rapidly, 
and several larg e efforts for pipeline processin g, source 
characterization (jPtak &: Griffithsll200a [ Aldcro ft 2006), 
and catalog generation are underway. The 2XMM cata- 
logue (Page 2006) is availabl e, and the Chandra so urce 
catalog". "'is due out in 2008 (jFabbiano et al.[[2007f) . 

Serendipitous wide-area surveys with Chandra were 
pioneered by the Chandra Multiwavelengt h Project 
(ChaMP; [Green et all [200l [Kim et al.[ |2004) for high 
Gala ctic latitude, while ChaMPlane (IGrindlav et al.l 
[2005( 1 has studied the stellar content of Galactic plane 
fields. The ChaMP, described in more detail below, also 
performs multiwavelength source matching and spectro- 
scopic characterization. These efforts are greatly aug- 
mented by lar ge surveys s uch as the Sloan Digital Sky 
Survey fSDSS: lYork et al.[[2000D . which will obtain spec - 
troscopyof- 100, 000 QSOs (e.g.. [Schneider et al.[[2007f l. 
and can additionally facilitate the efficient extrapolation 
of photometric quasar selection with photometric red- 
shift estimation al most a magnitude fa inter, towards a 
million QSOs (e.g.. [Richards et al.l[2008l ). 

The current paper studies the X-ray and optical prop- 
erties of the subset of these QSOs imaged in X-rays 
by Cha ndra as part of th e Chandra Multiwavelength 
Project (|Green et al.l[2004f ). 

2. THE QUASAR SAMPLE 

2.1. The SDSS Quasar Sample 

Most large samples of Type 1 QSOs are based on opti- 
cally selected quas a rs confirmed via optica l spectroscopy 
(|Bovle et al.l \198A [Schneider efall [1994 [Hewett et all 

^ F is the photon number index of an assumed power-law con- 
tinuum such that Ne{E) = . In terms of a spectral index 
a from = fugV°', we define F = (1 — a). 

^'^ OLox is the slope of a hypothetical power-law from 2500 A to 
2 keV; a^^ = 0.3838 log{l^^^^ -^/hkev) 

The Chandra source catalog 



|http : //cxc. harvard. edu/cscj 
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Il995l) . The largest, most uniform sample of optically se- 
lected quasars by far has been compiled from the SDSS. 
With the completion of the SDSS, we can expect some 
100,000 spectroscopically-confirmed quasars. The SDSS 
quasars were originally identified to i < 19.1 for spec- 
troscopy by their UV-excess colors, with later expansion 
for z > 3 quasars to i = 20.2 using ugri color crite- 
ria (|Richards et al.|[2"002[ ). The large catalog, the broad 
wavelength optical photometry, and subsequent foUowup 
in other wavebands has meant that research results from 
the SDSS spectroscopic quasar sample better character- 
ize the breadth of the quasar phenomenon than ever be- 
fore. However, the limited number of fibers available, 
fiber placement conflicts, and above all the bright mag- 
nitude limits of SDSS fiber spectroscopy mean that 10 
times as many quasars have been imaged, and could be 
efficiently identified from existing SDSS photometry. 

Using large spectroscopic AGN samples as 'training 
sets' can produce photometric classification and redshifts 
of far greater completeness and depth than spectroscopy. 
Without spectroscopic confirmation, photometric selec- 
tion criteria strike a quantifiable balance between com- 
pleteness and efficiency, i.e., a probability can be assigned 
both to the classification and the redshift. Deep pho- 
tometric redshift surveys like COMBO-17 (jWolf et al.l 
[2003) have found AGN to i? = 24 and z = 5 with high 
completeness. Efficient photometric selection of quasars 
in the SDSS using a nonparametric Bayesian classifi- 
ca tion based on kernel density estimation is described 
in [Richards et al.l (12004) for SDSS point sources with 
i < 21. An empirical algorithm to deter mine photometric 
redshifts for such quasars is described in lWeinstein et al] 
()2004[ ). The spectroscopic training samples for these 
methods now include far more high redshift quasars, and 
so the algorithms have been re-trained to include ob- 
jects redder than {u-g)=1.0 to classify high-z quasars, 
and applied to the much larger SDSS Data Release 6 
(DR6). This large catalog of ~1 million photometrically- 
identified QSOs and their photo metric redshifts is de- 
scribed in [Richards et al.l ()2008[ ). They estimate the 
overall efficiency of the catalog to be better than 72%, 
with subsamples (e.g.. X-ray detected objects) being as 
efficient as 97%. These estimates are based on an anal- 
ysis of the autoclu stering of the objects in the catalog 
(|Mvers et al.ll2006l ). which is very sensitive to stellar in- 
terlopers. However, at the faint limit of the catalog some 
additional galaxy contamination is expected. 

For luminosity and distance calculations, we adopt a 
Ho = 70 kms-^Mpc-i, Qa = 0.7, and = 0.3 
cosmology throughout. We assume a = —0.5 for the 
optical continuum power-law slope (v oc i^", where v 
is the emission frequency), and derive the rest-frame, 
monochromatic optical luminosity at 2500A (/ ° ; 
units erg s~^Hz~^) using the SDSS dereddened magni- 
tude with central wavelength closest to (1 + z) x 2500A. 

2.2. The Extended Chandra Multiwavelength Project 
(ChaMP) 

The Chandra Multiwavelength Project (GhaMP) 
is a wide-area serendipitous X-ray survey based on 
archival X-ray images of the {\b\ > 20deg) sky ob- 
served with the AXAF CCD Imaging Spectrometer 
(ACIS) onboard Chandra. The full 130-field Cy- 



cle 1-2 X-ray catalog is pubhc (|Kim et al.l l2007af) . 
and the most comprehensive X-ray number counts 
(logA^-logS) to date have been produced, thanks to 
6600 sourc e s and massive source-retrieval simulations 
(|Kim et al.l l2007bf ). We have also published early 
results of our deep (r ^25) NOAO /MO SAIC opti- 
cal imaging campaign (| Green eTall [2001 . now ex- 
tended to 67 fields (Barkhouse et al 2008, in prepa- 
ration). ChaMP results and data can be found at 
http : / /hea-www . harvard . edu/CHAMP, 

To improve statistics and encompass a wider range of 
source types, we have recently expanded our X-ray anal- 
ysis to cover a total of 392 fields through Chandra Cy- 
cle 6. We chose only fields overlapping with SDSS DR5 
imaging. To ease analysis and minimize bookkeeping 
problems, the new list of Chandra pointings (observation 
IDs; obsids hereafter) avoids any overlapping observa- 
tions by eliminating the obse rvation with the sh orter ex- 
posure time. As described in I Green et al.l ()2004f ). we also 
avoid fields with large (>,3') extended sources in either 
optical or X-rays (e.g., nearby galaxies MlOl, NGC4725, 
NGC4457, or clusters of galaxies MKW8, or Abell 1240). 
Spurious X-ray sources (due to e.g., hot pixel, bad bias, 
bright source readout streak s, etc.) have been flagged 
and removed as described in lKim et al. (2007a). Of the 
392 ChaMP obsids, 323 overlap the SDSS DR5 footprint. 

The GhaMP has also developed and implemented an 
xskycover pipeline which creates sensitivity maps for 
all ChaMP sky regions imaged by ACIS. This allows 
(1) identification of imaged-but-undetected objects (2) 
counts limits for 50% and 90% detection completeness 
and (3) corresponding fiux upper-limits at any sky posi- 
tion as well as (4) fiux sensitivity vs. sky coverage for any 
subset of obsids, as needed for logA^-logS' and luminosity 
function calculations. Our method is decr ibed in the Ap- 
pendix, and has been verified recently bv lAldcroft et aP 
fl200^) using the Chandra Deep Field South (CDF-S). 
The final sky ChaMP/SDSS coverage area (in deg^) for 
the 323 overlapping fields as a function of broad band 
{B band 0.5-8 keV) flux limit is shown in Figure [1] (see 
caption for specific definition of this limit). The area 
covered at the brightest fiuxes is 32deg.^ On average 5 
CCDs are activated per obsid. 

We have downloaded into the ChaMP database all the 
SDSS photometry, and the list of photo-z quasar can- 
didates within 20' of the Chandra aimpoint for each 
such obsid. Because the Chandra point spread func- 
tion (PSF) increases with off-axis angle (OAA), compar- 
atively few sources are detected beyond this radius, and 
source centroids also tend to be highly uncertain. Of X- 
ray detected candidates, we will show in Section l2T4l that 
98% of these candidates with spectra are indeed QSOs. 

Next we describe the identification of high-confidence 
ChaMP X-ray counterparts to SDSS QSOs in § O We 
then discuss in § 12.41 spectroscopic identifications for 
these objects. Section |3] then describes results for several 
interesting QSO subsamples, including our treatment of 
SDSS quasars that were not X-ray detected. 

2.3. X-ray/Optical Matching 

^■^ For 14 obsids, we extended to 28' radius, to achieve full cover- 
age of the Chandra footprint. For other obsids, the SDSS imaging 
strips do not completely cover the Chandra field of view. 
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Fig. 1.— Sky area vs B band (0.5-8 keV) flux limit for the 323 
obsids included in our ChaMP/SDSS field sample. Flux limit is 
defined here as the number of counts detectable in 90% of simula- 
tion trials, converted to flux assuming a power-law F = 1.7 at 2=0 
and the Galactic A'^f appropriate to each obsid. Chip S4 (CCD 
ID #8) is excluded throughout. The area covered at the brightest 
fluxes is 32 deg^ . 



The positional uncertainty of ChaMP X-ray source 
centroids has been c a refully analyzed via X-ray simu- 
lations by iKim et all (|2007b( ) and depends strongly on 
both the number of source counts and the OAA. Cross- 
correlating ChaMP X-ray source centroids with QSOs, 
we find that 95% of the X-ray/optical separations dxo of 
matched QSOs are smaller than the 95% X-ray positional 
uncertainty (Txp95- We first perform an automated 
matching procedure between each optical QSO position 
and the ChaMP 's X-ray source catalog. We adopt 4" 
as our matching radius criterion. Figure [5] shows that 
95% of the matched sample has an X-ray/optical posi- 
tion difference of <3", and expansion to include match 
radii above 4" would not substantially increase the sam- 
ple size. Searching the ChaMP catalog for X-ray sources 
within 4" of the optical SDSS quasar coordinate, yields 
1376 unique matches in the "Matched" sample. 

Although the overall efficiency of the SDSS pho- 
tometric QSO catalo g is only expected to be 72% 
(|Richards et al.l [2008D . the rarity of luminous Type 1 
quasars and X-ray sources means that matched objects 
should be quite clean. By repeatedly offsetting the SDSS 
coordinates of each QSO by 36" and rematching to the 
ChaMP X-ray catalog, we derive a spurious match rate 
of just 0.7%.^'^ This excellent result is thanks to Chan- 
dra's ^arcsec spatial resolution, which at SDSS depths 
allows for unambiguous counterpart identification, given 
that both Type 1 QSOs and X-ray sources are relatively 
sparse on the sky. 

We identify and remove a variety of objects with poten- 
tially poor data, including overlapping multiple sources 
(some of which are targeted lenses) , bright X-ray sources 
suffering from pile-up, optical sources with photometry 
contaminated by close brighter sources or within large 
extended galaxies or stellar diffraction spikes, etc. 

The ratio of the number of X-ray matches to the optical con- 
trol sample to the number of matches in the actual sample is 0.0067. 




Fig. 2. — Positional ofi'sets of matched sample. LEFT: The his- 
togram of X-ray/optical centroid separations in the full matched 
sample of 1376. The mode is 0.3", the median is 0.76", and the 
mean is 1.1". RIGHT: The ratio of the X-ray /optical centroid sep- 
aration to the 95% X-ray positional uncertainty vs. separation. 
The fraction of objects with separations larger than the 95% un- 
certainty (i.e., with ratio>l) remains relatively constant, and vir- 
tually no separations wider than 3" are larger than the positional 
uncertainty crxP95- 

In addition to the automated matching procedure, we 
also perform visual inspection (VI) of both X-ray and 
optical images, overplotting the centroids and their asso- 
ciated position errors. We retain only the highest confi- 
dence matches (matchconf =3). Most of the 105 objects 
we thereby eliminate have large ratios of dxo/<^xP95, or 
multiple candidate optical counterparts. We note that 
some of the most interesting celestial systems may be 
found among sources with matchconf <3. For example, 
these might include QSOs that are lensed, have bright 
jets, or are associated with host or foreground optical 
clusters or galaxies. Systems that are poorly-matched, 
multiply-matched, or photometrically contaminated may 
account for up to ~10% of the full X-ray-selected sam- 
ple. We therefore caution against blind cross-correlation 
of large source catalogs (e.g., the Chandra Source Cat- 
alog)^'* without such detailed quality control and visual 
examination of images. However, since we seek here to 
analyze the multiwavelength properties of a large clean 
sample of QSOs, and since most of these more compli- 
cated systems require significant further analysis or ob- 
servation, we defer their consideration to future studies. 

2.4. Spectroscopic Redshift Information 

After cross-correlation with the X-ray catalog, we 
sought spectroscopic redshifts for any objects in the pho- 
tometric QSO catalog. For this purpose, we obtained 
redshifts from existing ChaMP spectroscopy, from the 
SDSS (DR6) database itself, and then finally we searched 
the literature by cross-correlating optical positions with 
the NASA Extragalactic Database (NED), using a 2" 
match radius. Of 1376 matched objects, we found high 
confidence spectroscopic redshifts for 407, of which 43 
spectra were observed by the ChaMP. In striking tes- 
timony to the quality of the quasar selection algorithm 
(especially once candidates have been matched to X-ray 
sources) only 8 of these (2%) are not broad hue AGN. 
Three are narrow emission line galaxies, 2 are absorp- 
tion line galaxies (no emission lines of equivalent width 
Wx > 5A), and 3 are known BL Lac objects. We exclude 
these objects from all samples described below. Because 
the magnitude distribution (and therefore the photomet- 

^•^ The Chandra Source Catalog, available at 
'http : / /cxc . cf a. harvard. 6du/csc[ contains all X-ray sources 
detected by the Uhandra X-ray Observatory. 
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Fig. 3. — SDSS i magnitude vs. broadband (0.5-8 keV) X-ray flux 
for the SDSS/ChaMP QSO the Main sample. The flux shown is 
based on the best-flt PL-free yaxx model, with the effect of Galac- 
tic absorption removed from both X-ray and optical. X-ray de- 
tections are marked by black dots, and flux upper limits by green 
arrows. Open red triangles show QSOs with existing spectroscopic 
redshifts, clearly biased towards brighter optical mags. Radio-loud 
quasars (large open blue circles), BAL QSOs (open black squares), 
NAL QSOs (open black diamonds) and NLSls (open black trian- 
gles), and Chandra PI targets (large asterisks) are also indicated 
as shown in the legend. 

ric color errors) for objects lacking spectroscopic classifi- 
cations is fainter (see Fig [3]), we expect that the overall 
fraction of misclassified photo-z QSOs is larger than 2%, 
but it would be difficult to estimate it without deeper 
spectroscopic samples. A plot of the photometric vs. 
spectroscopic redshifts of these quasars is shown in FigO 

3. SAMPLES. SUBSAMPLES, AND DETECTION 
FRACTIONS 

The definition of our main sample and a variety of sub- 
samples are described in subsequent sections and sum- 
marized in Table [T] The size of our Main sample al- 
lows us to investigate the effects of luminosity or red- 
shift limits, X-ray non-detections, PI target bias, strong 
radio-related emission (RL QSOs), broad and narrow 
line absorption (BAL and NAL QSOs) , and Narrow-Line 
Seyfert Is (NLSls). In Table [1] samples to which we re- 
fer frequently are arrayed under 'Primary Samples' in 
decreasing order of the number of detections. Samples 
mentioned only once or twice in this paper are listed in 
similar order under 'Other Samples'. Tables listing bi- 
variate statistical results (Tables [5]- [7]) later in the paper 
list samples in this same order for reference. 

We define the Main sample to be the 2308 SDSS QSOs 
that fall on Chandra ACIS chips in a region of effective 
exposure > 1200 sec (excluding CCD #8; see below), re- 
gardless of Chandra X-ray detection status. Our cleaned 
matched sample (the MainDet sample) of X-ray/optical 
matched QSOs contains 1135 distinct X-ray sources with 
high optical counterpart match confidence, where we 
have removed all sources (1) with significant contami- 
nation by nearby bright optical sources, (2) significant 
overlap with other X-ray sources, (3) detected on ACIS- 
S chip S4 (CCD #8), because of its high background and 
streaking, (4) dithering across chips (which renders un- 
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Fig. 4. — LEFT: Histogram of effective Chandra exposure times 
for QSOs in the Main sample, in units of ksec, with bin size 6 ksec 
for detections (solid black line) and limits (dashed blue line). For 
detections, the mean and median exposure times are 25.9 and 
17.6 ksec, respectively. Inset shows detail at low exposure times 
using bin size 1 ksec. RIGHT: Histogram of Chandra off-axis an- 
gles for the Main sample, in units of arcmin. Black solid histogram 
shows how detections trend towards small OAA. Mean and median 
for detections are 6.4' and 5.8', respectively. Blue histogram shows 
that limits trend towards large OAA. Dashed histograms show the 
corresponding cumulative fractions. About 90% of the detections 
(compared to ~ 50% of the limits) are at OAA<12' off-axis. 

reliable the yaxx X-ray spectral fitting described below) 
or (5) with spectroscopy indicating that the object is not 
a Type 1 QSO. 

The mean exposure time for the MainDet sample is 
25.9 ksec per QSO, with an average of 3.6 QSOs de- 
tected on each Chandra field. For the 1173 non- 
detections in the Main sample, the mean exposure time 
is 17.6 ksec. A histogram of exposure times for all QSOs 
in the Main sample is shown in Fig[4l 

We publish key data for 1135 QSOs in the Main- 
Det sample in Tabled marking 82 sources that are the 
intended Chandra principal investigator (PI) targets. X- 
ray sources in plots include only the MainDet sample or 
subsets of it. Fig [6] shows luminosity vs. redshift for 
the MainDet sample. A large fraction of the z > 4 ob- 
jects are Chandra targets (large black stars). Strong 
redshift-luminosity trends are seen both in optical and 
X-ray, as is expected from any flux-limited survey. How- 
ever, the factor ~30-50 range in luminosity is unusual for 
a single sample; such breadth is usually only achieved 
using sample compilations encompassing diverse selec- 
tion techniques. In Fig [51 the large number of objects 
in our sample makes it difficult to distinguish the point- 
types presenting object class information, so Fig[7]shows 
a zoom-in on the most densely-populated regions of the 
L — z plane. 

Since we start with Type 1 SDSS QSOs, we are study- 
ing an optically selecte d sample, a.nd the selection func- 
tion is complex ( Richards et al.l l2006af) . If we limit 
the analysis to detections-only, then the sample is both 
optically- and X-ray-selected, and the selection function 
becomes increasingly complex. If instead we include 
all X-ray upper limits in the analysis, the sample re- 
mains fundamentally optically selected, but then statis- 
tical analyses must incorporate the non-detections (see 
§1111). 

The ChaMP 's xskycover pipeline allows us to investi- 
gate the detection fraction for the full SDSS QSO sam- 
ple, shown in Fig El Of 2308 SDSS QSOs that faU on an 
ACIS chip (the Main sample) in our 323 ChaMP fields. 

For 14 of 323 Chandra fields, no QSOs are detected. 
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1135 (49%) are detected in the MainDet sample. De- 
tection fractions as a function of SDSS QSO mag and 
redshift for this sample are shown in Fig [5) To minimize 
sample biases, we can also examine detection fractions 
as a function of X-ray observing parameters like expo- 
sure time and OAA. To simultaneously optimize detected 
sample size and detection fraction, we simply maximize 
N^^^/Niim, where N^et and Num are the number of X-ray 
detections and non-detections (flux upper limits), respec- 
tively. We find that an X-ray-unbiased subsample with 
a significantly higher detection rate is achieved by limit- 
ing consideration to the 1269 QSOs with OAA< 12' and 
exposure time T > 4ksec, of which 922 (72%) are detec- 
tions. This high detection fraction sample is called the 
D2L sample (Table [J). 

Data for X-ray non-detections is available in Table ??. 
We include in Table ?? only the 347 limits in the 
D2L sample, where the flux limits are sensitive enough 
to be interesting. By "interesting", we mean that the 
limits are close to or brighter than the faint envelope of 
detections. A flux limit several times brighter than that 
envelope provides no statistical contraints whatsoever on 
the derived distributions or regressions. For X-ray non- 
detections, data are more sparse all around for several 
reasons. QSOs with limits are optically fainter (mean 
and median i=20.4mag for the D2L sample limits com- 
pared to i = 19.9 for detections). Being fainter, fewer 
have SDSS spectroscopy. Also, as non-detections, none 
have been targeted for spectra by the ChaMP, so globally 
only about 10% of non-detected QSOs have optical spec- 
tra. The fraction of radio detections is also smaller (1.2% 
vs 4.8%). None are Chandra PI targets. Finally, X- 
ray non-detections lacking optical spectroscopy are some- 
what less likely to be QSOs. The selection efficiency 
(fraction of QSO candidates that are actual QSOs) be- 
tween about 0.8 < z < 2.4 is -95% (Richards ct al. 2004; 
iMvers et al.]l2006( l. but iRichards et all (j20d8h estimate 
that near the faint limit of i ~20.4mag, the overall QSO 
selection efficiency is ~ 80%. Particular attention must 
be paid to possible galaxy contamination at the faint 
end as the autoclustering estimates of the efficiency do 
not include galaxy interlopers at faint limits where SDSS 
star-galaxy separation begins to break down. However, 
many of these 'spurious' cross-matches may turn out to 
be (e.g., low-luminosity) AGN. In any case, the increased 
level of contamination by non-QSOs is another rationale 
for limiting the number of non-detections to those with 
sensitive X-ray limits. The nature of the statistical anal- 
ysis (discussed in § 16. ip is such that non-detections are 
included, but are assumed to follow the distribution of 
detections, and so effectively have a lower weight in the 
results. 

Amongst the undetected QSOs in the Main sample, 
165 have SDSS spectroscopy, of which 144 are high con- 
fidence spectroscopic QSOs. The 21 non-QSOs comprise 
16 stars, 5 galaxies. The higher (13%) rate of non-star 
spectroscopic classifications amongst undetected QSOs is 
not surprising, since X-ray detection greatly increases the 
probability that an optical AGN candidate is indeed an 
AGN. From the upper limit QSO sample, we remove the 
22 non-QSOs, and use the SDSS spectroscopic redshifts 
instead of the photometric redshifts wherever applicable. 





Fig. 5. — LEFT: Photometric redshift vs. spectroscopic redshift 
for QSOs detected in ChaMP fields. Filled (red) circles show QSOs 
for which the formal photometric redshift probability is >95%. A 
fraction of objects have large errors in their photometric redshifts. 
About 18% of QSOs with Zspec<l have Zphot > 2. This drops 
to 13% using only Pvoh{zphot)>0.5. RIGHT: Difference between 
the photometric and true (spectroscopic) redshift for QSOs in our 
sample, plotted against photometric redshift probability, again il- 
lustrates the reliability of these probabilities. 

3.1. Targets 

The "Non-target" detected sample (the noTDet sam- 
ple) of 1053 QSOs further eliminates 82 objects (7.2% 
of the MainDet sample) that are the intended targets of 
the Chandra observation wherein they are found. Tar- 
gets are on average brighter than most of the QSO sam- 
ple (see Fig |3]), but more importantly were chosen for 
observation for a variety of reasons unrelated to this 
study. In particular, targets tend to be more lumi- 
nous than serendipitous QSOs (Figs [6l [7]), and several 
are known lenses (e.g., HS 0818-f 1227, PG 1115+080, 
UM425=QSO 1120+019). The bias in sample character- 
istics is largely mitigated in the subsamples excluding 
targets (see Table [T|). However, the exclusion of tar- 
gets also produces a (much smaller) bias because some 
objects with similar characteristics would have been in- 
cluded (at a lower rate) were the Chandra pointings all 
truly random. Because many of the targets are indeed 
of interest (e.g., high-z QSOs), we include them in most 
discussions, but always check that results are consistent 
without them. We also note that some target bias proba- 
bly affects the X-ray sample even after the exclusion of PI 
target QSOs, because PI targets may cluster with other 
categories of X-ray sources such as other AGN, galaxies 
or clusters. Overall, a comparison of regression results^^ 
for several of our subsamples that differ only in target 
exclusion does not indicate a significant target bias, due 
at least in part to our large sample sizes. 

3.2. Radio Loudness 

Quasars with strong rad io emission are observed to be 
more X-ray luminous fe.g.. lGreen et al.lll995l : [SIien et al.l 
12006( 1. At least some of the additional X-ray luminosity 
is likely to originate in physical processes related to the 
radio jet rather than to the accretion disk, so it may be 
important to recognize those objects that are particularly 
radio loud. 

The Faint Images of the Radio Sky at Tw enty- 
Centimeters (FIRST) survey (jBecker et al.l |1995[ ) from 
the NRAO Very Large Array (VLA) has a typical (5(t) 
sensitivity of —1 mJy, and covers most of the SDSS foot- 
print on the sky. Following llvezic et al.l (|2002[ ). we adopt 

Table l6l7l shows similar results comparing e.g., the Main- 
Det sample and the noTDet sample, or the D2LNoRB sample and 
the D2LNoTRB sample. 
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Fig. 6.— Luminosity at 2500A (LEFT) and 2keV (RIGHT) vs. redshift for the SDSS/ChaMP the Main sample. Above z ~ 2.5, the 
number of QSOs dechnes steeply, due to the SDSS ma gnitude l imit and to the decreased efficiency of the photometric selection algorithm 
as it crosses the stellar color locus (see Fig |8] and Richards et al.i 120021 ). These plots show 56 QSOs with 2 > 3, of which 34 are new 
serendipitous detections. X-ray upper limits are shown as small green triangles here. See Fig [3] for symbol types. The dashed black 
rectangle surrounds the zBoxDet sample, a portion of the '2500 A ~^ plane chosen to test for redshift dependence. The green rectangle 

"LoptBox" surrounds the LoBox sample, used to test for dependence on Z » . The blue rectangle surrounds the zLxBox sample, chosen 

2500 JT- 

to avoid X-ray flux limit bias. 
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Fig. 7. — LEFT: Zoom-in of 2500A luminosity vs. redshift. Objects with spectroscopic redshifts (open red triangles) tend to be at high 
optical luminosities by selection. Detectably radio-loud QSOs are shown with open blue circles. See Fig |3] for symbol types. BALQSOs 
(open black squares) are mostly detected at 2: > 1.6 where the CIV region enters the optical bandpass. There appears to be no preference 
of BALQSOs for high optical luminosity, apart from the bias caused by Chandra target selection. RIGHT: Zoom-in of 2keV luminosity 
vs. redshift. Here the RL QSOs clearly populate the upper luminosity envelope. BALQSOs are preferentially X-ray quiet, unlike the QSOs 
with NALs only (open diamonds). 



a positional matching radius of 1.5", which should result 
in about 85% completeness for core-dominated sources, 
with a contamination of ~3%. We thereby ma tch 69 
sources to the Matched sample. iJiang et al.i ()2007l ) 
matched the FIRST to SDSS spectroscopic quasars and 
found that about 6% matched within 5". We might ex- 
pect a lower matched fraction because our optical pho- 
tometric sample extends 1-2 mag fainter. On the other 
hand, we are looking at X-ray-detected quasars, so the 
actual matched fraction of ~ 5% is similar. We also 



matched all the quasar optical positions to the FIRST 
within 30", and visually examined all the FIRST im- 
ages to look for multiple matches and/or lobe-dominated 
quasars. There are 26 sources we judged to have reliable 
morphological complexity that are resolved by FIRST 
into mult iple sources. The NRAO VLA Sky Survey 
(NVSS; iCondon et al.lll998f ). lists detections for 19 of 
these. Comparing NVSS to summed FIRST fluxes, we 
found NVSS fluxes slightly larger - < 2% difference in 
the mean 20% max). Since the NVSS beam is larger 
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Fig. 8. — LEFT: Histogram of redshifts for detected (solid blue), 
non-detected (red dashed) and all (black solid) QSOs in the the 
MainDet sample. We tally the 'best' redshift for each object (i.e., 
spectroscopic redshifts are always used when available). Above 
z ~ 2.5, the number of QSOs declines steeply, due to the decreased 
efficiency of the photometric selection algorithm as it crosses the 
stellar color locus (Richards et al. 2002). The inset shows the de- 
tected fraction as a function of redshift. RIGHT: Histogram of 
SDSS i mag for detected (solid blue), non-detected (red dashed) 
and all (black solid) QSOs. The inset shows the detected fraction 
as a function of magnitude. 
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Fig. 9. — LEFT: Histogram of redshifts for detected (solid blue), 
non-detected (red dashed) and all (black solid) QSOs, after re- 
striction to obsids with T > 4ksec and QSOs with OAA<12' (the 
D2L sample). The inset shows the detected fraction as a function 
of redshift. RIGHT: Histogram of SDSS i mag for detected (solid 
blue), non-detected (red dashed) and all (black solid) QSOs. The 
inset shows the detected fraction as a function of magnitude. 



(45") than FIRST (5"), FIRST detection algorithms may 
exclude some of the extended source flux as background, 
so we include the NVSS fluxes for these 19 objects, and 
summed FIRS T fluxes for thos e remaining. 

Following I vezic et al.l 1)200 2'). we adopt a radio loud- 
ness parameter R as the logarithm of the ratio of the ra- 
dio to optical monochromatic flux: R = \og{F20cm/Pi) = 
OA(i — 17120 cm), whe re TO2o cm is the radio AB magnitude 
(|Oke fc Gunnl lM. msocm = -2.5 log(F2o c™/3631 Jy) 
calculated from the integrated radio flux density, and i is 
the SDSS i-band magnitude, corrected for Galactic ex- 
tinction. We adopt a radio-loudness threshold i? = 1.6. 
Thus there are 72 QSOs in the Main sample with radio 
detections, of which 57 (79%) are radio-loud. For the 
MainDet sample (detections only), there are 55 radio- 
detected QSOs, of which 43 (78%) are radio-loud. 

Given the ('^l mJy) source detection limit of the 
FIRST Survey, all RL QSOs will be detected to about 
i ~20.4mag. For the magnitude range 17< i <20 where 
the statistics are good and the FIRST is sensitive to all 
RL QSOs, we find 41 of 529 (8±1%) such QSOs from 
the MainDet sample are detected by the FIRST, with 29 
(5.4%) that are radio loud. Since the 34% of our full the 
MainDet sample that is fainter than j=20.4 suffers from 
incomplete radio loudness measurements, some 2% may 
be unidentified RL QSOs. A similar fraction pertains if 
we count X-ray non-dctcctions as well (the Main sample) . 




1 2 

Redshift 

Fig. 10. — Radio-loudness vs. redshift for the MainDet sample 
(detections). Radio-loud objects R > 1.6 are shown with open blue 
circles. Radio-quiet but FIRST radio-detected objects are shown 
as filled green circles. All other symbols (described in FigO have 
radio flux upper limits only. Note that most Chandra targets are 
distinctly either loud or quiet, highlighting a bias in the target 
subsamplcs. 

3.3. Broad and Narrow Absorption Line Quasars 

We identified QSOs with broad absorption lines 
(BALs) and narrow absorption lines (NALs) directly by 
visual inspection of QSOs with spectroscopy, finding 16 
BAL and 11 NAL QSOs in the the MainDet sample. Ten 
(two) of the BAL (NAL) QSOs were the Chandra PI tar- 
gets. 

The best estimates to date of the raw BALQSO frac- 
tion amon g optically selected quasars range from about 
13 - 20% (jReichard et "311120031 : iHewett fc Foltg r2003). 
From the SPSS D R3 sample of iTrump et"aI1 p006) . 
iKnigge et all (|2008[ ) carefully define what is a BALQSO 
and correct for a variety of selection effects to derive an 
estimate of the intrinsic BALQSO fraction of 17%±3. 
The vast majority of BALQSOs in the SDSS are above 
redshift 1.6 because only then does the CIV absorption 
enter the spectroscopic bandpass. If we determine our 
BALQSO fraction in the MainDet sample only counting 
the serendipitous (non- Target) QSOs with z > 1.6 and 
spectroscopic redshifts, we find just 4 out of 119 QSOs 
with BALs. Even with sensitive X-ray observations such 
as these, X-ray selection is strongly biased against the 
high ionized absorbing columns along the line of sight 
towards the X-ray emitting regions of BALQSOs. Of the 
24 absorbed (BAL or NAL) QSOs two are detectably ra- 
dio loud; SDSS J171419. 24-1-611944.5 - a BAL QSO - and 
SDSS J171535.96-f632336.0 - a NAL QSO - are targets 
selected (by Chandra PI Richards) as reddened QSOs. 

3.4. Narrow Line Seyfert Is 

X-rays from NLSls are of particular interest because 
they were thought to sho w marked variability and strong 
soft X -ray excesses (e.g., ([Green et al.lll993HBoller et alj 
Il996l )). NLSls are proposed to be at one extreme of the 

A much smaller number of the rare low-ionization BALQSOs 
(with BALs just blueward of Mgll) are found at lower redshifts. 
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so-called (jBoroson fc Greenlll992D "eigenvector 1" , which 
have been suggeste d to correspond to low SMBH masses 
(|Grupe et a l."20d?) and/or high (near-Eddington) accre- 
tion rates (Boroson 2002). 

For objects with SDSS spectra encompassing H/3 
z < 0.9), we identify as Narrow Line Seyfert Is 
Osterbrock fc Poggdri985t NLSls hereafter) those ob- 
jects with FWHM(H/3)<2000 km s'^ and line flux ratio 
[0III]/H/3< 3. The FWHM measurements are obtained 
via FWHM=2.35 c(t/Ao/(1-I-2:), where is the variance 
of the Gaussian curve that fits the H/3 emission line, z is 
the redshift, and Aq is the rest- frame wavelength of the 
H/3 line (4863 A) . We extract measurements of a (in A) 
and line fluxes from the SDSS DR6 SpecLine table. 

Our MainDet sample contains at least 19 NLSls (those 
with spectroscopy in the redshift range to include H/3). 
Among the non-detections, we identify 3 more. Six of 
these 22 objects are targets described in (jWilliams et al.l 
120041 ) ■ and SDSS J125140.33-t-000210.8 was a target se- 
lected (by Chan dra PI R ichards) as a dust-reddened 
QSO. The Williams et al.l (|2004 ) sample consisted of 17 
SDSS NLSls selected from the RASS to be X-ray-weak. 
Their study confirmed earlier suggestions that strong, ul- 
trasoft X-ray emission is not a universal characteristic of 
NLSls. 

We therefore present here new results for a sample of 
15 SDSS NLSls observed by Chandra. This admittedly 
small sample is nevertheless the largest published sample 
of optically-selected NLSls with unbiased X -ray observa- 
tions. (The sample of I Willia ms ct al. (2004) was selected 
to be X-ray weak, while the Grupe et al. (2004) ROSAT 
sample was selected to have strong soft X-ray emission.) 
We find no evidence for unusual SEDs from the distribu- 
tions of either aox or F. 

4. OPTICAL COLORS AND REDDENING 

In Fig [Til we plot the {g — i) colors of the matched 
SDSS/ChaMP QSO sample as a function of redshift, and 
compare to the optical-only sample.^'*. The Chandra- 
detected sample does not show significantly different col- 
ors from the full optical sample. This likely attests to 
(1) the sensitivity of the Chandra imaging relative to 
the magnitude limit of the optical sample and (2) the 
fact that Type 1 QSOs are largely unabsorbed in both 
the optical and X-ray regimes. 

The right panel of Fig [Tl] shows that 10 of 14 
BALQSOs are above the A(g-i)=0 line. This reflects 
that SDSS BALQSOs tend to be redd er than average 
(jReichard et al.l l2003t iDai et all l2007f l. Most of the 
RL Q SOs are also redder than average. 'Richard s et al.l 
()2001l ) found a higher fraction of intrinsically red- 
dened quas a rs arn ong those with FIRST detections. 
Ilvezic et al.l (|2002D found that RL QSOs are redder 
than the mean (at any given redshift) in (g — i) by 
0.09±0.02mag. Fig [Til confirms a similar trend in the X- 
ray detected SDSS/ChaMP sample. At the same time, 
a small number of RL QSOs are found on the blue ex- 
treme of the color-excess distribution. These trends are 

18 The optical-only sample refers to all SDSS QSOs within ~20' 
of all Chandra pointings, regardless of whether its position falls on 
an ACIS CCD. We use a 9th-order polynomial fit to the optical- 
only sample with the following coefficients: 0.698892, 3.011733, - 
20.358267, 37.850353, -33.617121, 16.652032, -4.861889, 0.833027,- 
0.077526, 0.003026 



fully consistent with the detailed res ults found by stack - 
ing FIRST images of SDSS quasars ([White et al.l[200l . 
independent of X-ray properties. 

5. X-RAY SPECTRAL FITTING WITH YAXX 

Besides comparing the broadband multiwavelength 
properties of QSOs, Chandra imaging provides X-ray 
spectral resolution capable of yielding significant con- 
straints on the properties of emission arising nearest the 
supermassive black hole. While the ChaMP calculates 
Hardness Ratios (HR) and appropriate errors for every 
source, these can be difficult to interpret, since HR con- 
volves the intrinsic quasar SED with telescope and in- 
strument response, and does not take redshift or Galac- 
tic column N^"'^ into account. A direct spectral fit of the 
counts distribution using the full instrument calibration 
and known redshift and /V|^"' provides a much more di- 
rect measurement of quasar properties. Note that even 
in the low-counts regime one can obtain robust estima tes 
of fit parameter uncertainties using the Iclashl (fl979h fit 
statistic. 

We use an automated procedure to extract the spec- 
trum and fit up to three models to the data. For all 
objects in the Matched sample, we first define a circu- 
lar source region centered on the X-ray source which 
contains 95% of 1.5 keV photons at the given off-axis 
angle. An annular background region is also centered 
on the source with a width of 20". We exclude any 
nearby sources from both the source and background re- 
gions. We then use CIAO^^ tool psextract to create 
a PH A (pulse height amplitude) spectrum covering the 
energy range 0.4-8 keV. 

Spectral fitting is done using the CIAO Sherpc?^ tool in 
an automated script known as yaxx^^ (Aldcroft , ,2006). 
All of the spectral models contain an appropriate Galac- 
tic neutral absorber. For all sources we first fit two 
power-law models which include a Galactic absorption 
component frozen at the 21 cm value:^^ (1) fltting photon 
index F, with no intrinsic absorption component (model 
"PL") and (2) fitting an intrinsic absorber with neutral 
column A^IP*'" at the source redshift, with photon index 
frozen at F = 1.9 (model "PLfix"). Allowed fit ranges 
are -1.5 < F < 3.5 for PL and lO^^ < A^lf*'' < lO^^ for 
PLfix. These fits use the Powell optimization method, 
and provide a robust and reliable one-parameter charac- 
terization of the spectral shape for any number of counts. 
Spectra with less than 100 net counts^ "^ were fit u sing the 
ungrouped data with Cash statistics ([Cash[ll9"79h . Spec- 
tra with more than 100 counts were grouped to a mini- 
mum of 16 counts per bin and fit using the statistic 
with variance computed from the data. 

Finally, X-ray spectra with over 200 counts were also 
fit with a two-parameter absorbed power-law where both 
F and the A^^"' were free to vary within the above ranges 
(model "PL_abs"). 

1^ http://cxc.harvard.edu/ciao 

■^^ http://cxc.harvard.edu/sherpa 

http:/ /cxc. harvard.edu/contrib/yaxx 

■^^ Neutral Galactic column density A'^jj"' taken from 
IDickev &: Lockman I II1990I ) for the Chandra aimpoint position on 
the sky. 

Source counts derived from yaxx may differ at the ~1% level 
from those derived by ChaMP XPIPE photometry, due to slightly 
different background region conventions. 
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Fig. 11.— LEFT: SDSS (g - i) color vs. best redshift. Symbols show individual QSOs in the SDSS/ChaMP the MainDet sample. The 
black line shows the mean color at each redshift (in redshift bins of 0.12 for z <2.5 and 0.25 for higher redshifts) for the full optical QSO 
sample (regardless of Chandra imaging). We derive a smooth (9th order) polynomial fit to those means, whose value is plotted as the 
"expected" (g — i) with a green asterisk (at the redshift of each actual QSO). RIGHT: The color excess A(g — i) (the difference between 
the actual and "expected" (g — i) color) is plotted against redshift for a limited redshift range, to show highlights. The residuals of the 
polynomial fit to the mean binned (g — i) of the full optical sample are shown connected by a solid green line. See Fig|3]for symbol types. 
Most targets, BALQSOs, and RL QSOs are redder than average. 



5.1. X-ray Spectral Continuum Measurements 

We compile "best-PL" measurements, where for fewer 
than 200 counts, we use F from the PL {N^^*-"^ fixed at 
zero) fits and for higher count sources we use PL_abs 
(both r and A^IP**" free). In the MainDet sample there 
are 156 sources with 200 counts or more. High-count ob- 
jects are found scattered at aU luminosities below z '^2.5. 
QSOs with >200 counts (0.5-8 keV) with both T and 
]\[mtr g.j-g jj-^ yaxx, are weh-distributed in I2 keV amongst 
the detections, due to the wide range of exposure times. 

The mean F for all the 1135 QSOs in the MainDet sam- 
ple is 1.94±0.02 with median 1.93. Means and medi- 
ans for the MainDet sample and the subsamples dis- 
cussed in this section are listed in Table IH The typ- 
ical (median) error in F, AF ^ 0.5, is similar to the 
dispersion 0.54 in best-fit values of F. If we limit the 
sample to the 314 sources with >100 counts, the typi- 
cal error is 0.32, with no change in mean or median F. 
The F distribution we find is similar to that found re- 
cently for smaller samples of BLAGN. Just et al. (2007|) 
studied a sample of luminous optically-selected quasars 
observed by Chandra, ROSAT, and XM M-Newton, and 
found <F> = 1.92±0.09 for 42 QSOs. (|Mainieri et all 
|2007() studied a sample of 58 X-ray selected BLAGN in 
the XMM-COSMOS fields and fo und <F>=2.09 with a 
dispersion of ~0.26. iPaed (|2006D found <F>=2.0±0.1 
with a dispersion of ~0.36 for 50 X-ray selected BLAGN 
in the 13^ XMM-Newton/ Chandra deep field. 

Figure [12] shows a histogram of best-fit power-law 
slopes for several interesting subsamples of QSOs from 
both the MainDet sample and from the noTDet sample 
which omits PI targets. The mean and median values 
for these subsamples are listed in Table [H We do not 
separately plot the radio-quiet (RQ) QSO sample, since 
it follows quite closely the shape of the full sample his- 
togram. 



For 43 detected RLQSOs in the MainDet sample, 
the nominal mean slope is <F> = 1.73±0.05 with me- 
dian 1.65, with a distribution significantly flatter than 
for the 704 definitively RQ QSOs (<F> = 1.91±0.02) in 
the MainDet sample, using the two-sample tests de- 
scribed in § 16.1.11 below. RL QSOs are known to have 
fiatter high energy con tinuua from previous work (e.g., 
iReeves fc Turnerll2000[ ). 

For the 15 known BALQSOs in the MainDet sample, 
the nominal mean slope is <F>=1.35±0.15 with median 
1.3, and the distribution is significantly different than 
for the full MainDet sample (minus BALs, NALs, and 
NLSl), using the two-sample tests described in § 16.1.11 
below. The difference is even more significant {Pmax < 
0.01%) when comparing only to the 667 definitively RQ 
QSOs {RL < 1.6). Since the mean (median) number 
of (0.5-8 keV) counts for BALQSOs is just 27 (15), the 
nominally lower F likely reflects undetec ted absorption 
(|Green et al.l[2?)?nl : IcTallagher et al.l[200l ). 

The mean slope <F>=1.68±0.13 with median 1.75 for 
the 9 known NALs in the MainDet sample is indistin- 
guishable from the full MainDet sample (minus BALs, 
NALs, and NLSl), probably but the NAL statistics are 
poor. On the other hand, a smaller, non-overlapping 
sample of NALs o bserved by Chandra published by 
iMisawa et all (|2008f ) agrees that the X-ray properties of 
intrinsic NAL quasars are indistinguishable from those 
of the larger quasar population. 

For the 19 known NLSls in the MainDet sample, 
the nominal mean slope is <F>~2.01±0.15 with me- 
dian 1.95, indistinguishable from the comparison sample 
(MainDet sample minus BALs, NALs, and NLSl). 

5.1.1. X-ray Spectral Evolution 

No signs of evolution have been detected for the 
intrinsic power-law slope F of QSOs: z>A samples 
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Fig. 12. — LEFT: Histogram of best-fit X-ray spectral power-law slope F (best-PL) for the MainDet sample. The full sample (black solid) 
histogram has been divided by 15, and the RL QSO (long-dashed black) histogram by 2, for ease of comparison with smaller subsamples. 
BALQSOs (blue down-slash shading) show very flat slopes, due to strong intrinsic absorption. The NAL distribution (green up-slash 
shading) suggests a possible bimodality. Neither the NAL nor the NLSl (magenta dense shading) nor the hig redshift (short-dashed red) 
histograms are signiflcantly different by K-S test from the full sample distribution. RIGHT: Same plot, but for the noTDet sample which 
omits PI targets. The RL QSO distribution is less distinct here, and the softest (largest F) NLSls disappear. 

There is conflicting eviden ce as to whether F correlates 
with X-ray luminosity Lx (|Dai et al.l l2004l: iPage et all 
'2004), or pe rhaps w i th Eddington ratio fPorquet et al.l 
■ 2004 - Kellvl l2007hl : iShemmer et all [2008 : Kellv^eTan 
I2OO8I ). ISaez ()2008( ) find significant softening (increase) 
of F with luminosity in selected redshift ranges, for a 



(|Vignali et al.l [2005 IShemmer et"all [20061) show r^2. 
just like those at lower redshifts ([Reeves fc Turnei|[2000l ) . 
In a recen t small sample of high (optical) luminos- 
ity QSOs ([Just et al.l l2007l ) also found no trend of F 
with redshift. A larger compilation also shows at best 
marginal signs of evolution, and only for well-chosen red- 
shift ranges (Sacz 2008). 

Fig [T3| shows F vs. redshift both for the individual 
QSOs (top), and for subsamples (bottom) in fixed bins 
of A 2;. We also try another binning, where we grow 
redshift bins from z = until each bin has between 100 
and 150 objects (except for the high redshift bin z >3 
which has just 58). The largest difference between any 
of the bins is 0.17, whereas the typical rms dispersion in 
all bins is '^0.5. 

The noTDet sample and the NoRBDet sample remove 
contributions from targets, and from RL or absorbed 
QSOs, respectively, the HiCt sample is the subset of 
the MainDet sample with counts>200, where fits for 
N^'-*'^ are also performed. This decreases the effect of 
QSOs that may appear to have hard F but are actu- 
ally absorbed, the HiCtNoTRB sample also has only 
counts>200, but further removes targets, RL, BAL, and 
NAL QSOs. We find no evidence for evolution; the 
null hypothesis (no correlation) cannot be rejected with 
P < 8% for any of samples 0, 1, 2, 20, or 21, which 
test possible contamination or biases from targets, RL, 
and/or absorbed QSOs. 

Given that the typical uncertainty in F for a single 
QSO (~0.5) is comparable to the sample dispersion, si- 
multaneous spectral fitting o f subsamples of d etected 
QSOs in bins of redshift (ala iGreen et al.l l2001l ) might 
improve constraints, but such a project is both beyond 
the scope of this paper and probably unwarranted given 
the scarce evidence for evolution to date. 



5.1.2. X-ray Spectral Trend with X-ray Luminosity 



combination of Type 1 and Type 2 X-ray selected AGN 
from the CDFs. Type 2 objects dominate numerically, 
especially at low Lx, and show larger scatter in F at all 
luminosities, so Type 2 objects dominate the trends that 
they observe. Our sample is complementary in that it 
treats only Type 1 QSOs, and provides the largest, most 
uniform such spectral study to date. 

We detect a significant anticorrelation between F and 
hkev across a several subsamples (see Table [5] and 
Fig [Ml) , but none between F and either redshift or 



/ ° . The best-fit continuum slope F hardens (de- 
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creases) with increasing luminosity for the MainDet sam- 
ple. Since F is itself used to calculate I2 kev via the stan- 
dard X-ray power-law fc-correction, there might be some 
danger that the apparent anticorrelation is induced. We 
investigate this in two ways. First we examine the range 
in the ratio of I2 keV calculated assuming our best-fit F fc- 
correction to that calculated with a fixed T — 1.9. Across 
the full X-ray luminosity range, the different assump- 
tions induce a change of ~0.1dex (~ 25%), insufficient 
to account for the observed trend. Second, we examine 
whether F correlates significantly with I - , which is 

2500 A 

clearly unaffected by the X-ray fc-correction. There is 
indeed a significant trend in the MainDet sample (see 
Table [5]). 

One possibility is that the observed trend of F with 
luminosity might be due to an undetected increase in 
^mtr^ To test this, we examined the relationship in the 
HiCtNoTRB sample, the subset of the MainDet sample 
with counts>200, where fits for N"^*^^ are also performed. 
We do not claim that all absorption is detected in the 
HiCtNoTRB sample. However, if indeed undetected ab- 
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sorption caused the trend in the MainDet sample, we 
would expect the anticorrelation to weaken or flatten in 
the HiCtNoTRB sample. Instead, the anticorrelation is 
still significant, and the slope is actually steeper (see Ta- 
ble Ej). 

It is also conceivable that soft X-ray emission associ- 
ated with starformation activity might contaminate the 
sample at low luminosities. To test this, we create the 
HiLoLx sample with definitively QSO-like luminosities in 
both wavebands, i.e., log / °> 29.8 and hkeV > 26. 

' ' ° 2500 A ^ ^ ^ 

For this sample, the anticorrelation remains strong and 
and steep (see Table [5]) . 

We speculate that the balance between thermal (accre- 
tion disk) and non-thermal X-ray emission may shift to- 
wards higher luminosities. If at high Z2 fceV J non-thermal 
emission represents a larger fraction of the emitted X- 
rays, we would expect a hardening (decrease) of F with 
luminosity, as is seen. We would also expect to find more 
RL QSOs at high hkeV: and we do. Since quasars such 
as RL QSOs with a larger fraction of non-thermal emis- 
sion also show stronger X-ray emission relative to op- 
tical (smaller aox), we might expect a correlation be- 
tween F and aox such that as the spectrum hardens, 
aox decreases. This is shown to be the case in § 17.21 
below. The idea of a significant non-thermal, probably 
jet-rel ated emission comp o nent e ven in RQ QSOs is not 
new. iBlundell fc BeaslevI ()1998l ) found strong evidence 
from Very Long Baseline Array (VLBA) observations for 
jet-producing ce ntral engines in 8 of t he 12 RQ QSOs 
in their sample. iBarvainis et al.l (|2005| ) find similarities 
between RL and RQ quasars spanning a host of radio 
indicators: variability, radio spectral index, and VLBI- 
scale components, suggesting that the physics of radio 
emission in the inner regions of all quasars is essentially 
the same, involving a com pact, partially opaque core to- 
gether with a beamed iet. ICzernv et al] ([2008) similarly 
find evidence for a blazar component in RQ QSOs by 
modeling their variability. 

More sensitive empirical tests of whether the observed 
trend is due to substantial (but not directly detectable) 
absorption depressing the observed soft X-ray continuum 
at high luminosities, or to an increasing thermal frac- 
tion at lower luminosities could be performed by stacking 
counts in narrow energy bands from Chandra images of 
all QSOs (detected or not). Stacking h as been used ef- 
fectively this way in the CDFs fe.g.. ,Steffen et al.ll2007l : 
iLehmer et al.l 120071 ) . but the task is more daunting for 
the ChaMP, where care must be taken to properly ac- 
count for the effects of a much wider range of7V|"', CCD 
quantum efficiency values, and exposure tim es. 

A study of 35 Type 1 QSOs bv .Shemmer et al.l (|2008l ) 
finds that F increases (softens) with L/LEdd, the latter 
derived from FWHM(H/3) and i^L„f5100A). iKellv et all 
(|2008h examine a larger sample, and find complicated 
relationships between Mbh , Ljjv / LEdd, and Lx/LEdd 
that change direction depending on the emission line used 
to estimate Mbh ■ In a subsequent paper, we are extend- 
ing our current analysis to include estimates for Mbh 
and L/LEdd for the spectroscopic subsample of the cur- 
rent paper. 



5.2. X-ray Intrinsic Absorption Measurements 
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Fig. 13. — TOP: Best-fit X-ray spectral power-law slope F vs. 
redshift for tiie MainDet sample. Large open green circles show 
QSOs with more than 200 (0.5-8 keV) counts and simultaneous 
r-Nff*'~ fits. No evidence is seen for evolution in F. See Fig \3\ 
for symbol types. Dots to distinguish detections from limits in 
other plots are omitted here, since by definition all objects are 
detected. Radio-loud QSOs (open blue circles) tend to have fiatter 
slopes. Some of the fiattest slopes seen are for BALQSOs (open 
black squares), due to their strong intrinsic absorption. Typical 
mean errors on F are ~ 0.5 below 2 = 2, rising to 0.8 at the 
highest redshifts. BOTTOM: Mean F values are shown at the 
mean redshift for QSOs in 5 redshift bins of width A 2 = 1 for the 
MainDet sample. Errorbars show the error in the means for both 
axes. 
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Fig. 14. — TOP: Best-fit X-ray spectral power-law slope F vs. 
X-ray luminosity for the MainDet sample. See Fig |3] for symbol 
types. Large open green circles show QSOs with more than 200 
(0.5-8 keV) counts and simultaneous F-A''|p''" fits. A significant but 
shallow trend towards harder spectra (smaller F) is evident from 
the best-fit regression lines shown: the MainDet sample (red solid), 
the HiCtNoTRB sample with counts>200 (blue short-dashed) or 
HiLoLx sample with log I °> 29.8 and l2keV > 26 (blue long- 

dashed). Fit parameters are shown in Table [5] BOTTOM: Mean 
F values are shown at the mean luminosity for QSOs in 6 bins of 
width Alx = 0.75 for the MainDet sample. Errorbars show the 
error in the means for both axes. 
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The quaUty of X-ray measurements of intrinsic absorb- 
ing columns depends strongly on the number of counts 
in the quasar, but also on the redshift of the object, as 
illustrated in a simple test in FigfTSl 

Fig [16] shows our best-fit intrinsic absorption column 
measurements, which are overwhelmingly upper limits, 
either when assuming (as for counts<200) that T = 1.9 
or when fitting both T and N*^**^ (for QSOs with >200 
counts). Just one of the RL QSOs has detectable N^**", 
which is nevertheless not large (10^^cm~^). Jet-related 
X-ray emission may reduce any absorption signatures. 

Strong intrinsic absorption is relatively rare in radio- 
quiet Type 1 QSOs. In the strictest interpretation 
of AGN unification models, none of these broad line 
AGN should be significantly X-ray-absorbed. Small ob- 
scured fractions might be expected by selection, which 
requires that the broad line region (BLR) not be heav- 
ily dust-reddened, i.e. our view of the BLR is unob- 
scured. X-ray absorbed BLAGN have therefore some- 
times been called "anomalous" . The obscured fraction of 
BLAGN from the literature spans a wide variety of sam- 
ples and analysis methods, but most define "obscured" as 
N^f>iOffcnr:fand find fractions o f about 1 0% or less 
(jPerola et al.ll2004l : Page 2006; Maini eri et al . 2007). By 
contrast, selecting optically unobscured AGN only from 
optical/IR p hotometry, a nd us ing X-ray hardness ratios 
from XMM, iTaier et all (|2007l) find 30%. 

To estimate the obscured fraction of SDSS QSOs from 
the ChaMP, we first limit to the D2L sample and z < 2 
to maximize the fraction that are X-ray detected (see 
Fig [9]). We further limit to QSO-like optical luminosi- 
ties log / ° > 29.8, yieleding 630 QSOs, for which 498 
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(79%) are X-ray detections. If we rather stringently re- 
quire that the 90% lower bound to the best-fit log iVIJ'*'' 
exceeds 22, then 50 of the 498 (10%) of X-ray detected 
QSOs are obscured. If instead we simply use the best-fit 
7Vi^'*^ 98 of 498 (20%) of detected QSOs are obscured. 
However, if we make the unlikely assumption that all X- 
ray-undetected QSOs are missed due to heavy intrinsic 
absorption, the obscured fraction could be as high as 29% 
and 36% for these two different methods, respectively. 

We note that the correlation of dust to total column 
of X-ray absorbing metals is not strict. For instance, 
we know that a si gnificant fraction of optically selected 
QSOs ri3 - 20%: iReichard et all [20031: iHewett fc Foltd 
I2003t iKnigge et all l2008f) appear as BALQSOs, which 
are thought to be seen along a sightline piercing warm 
(ionized) absorbers. BALQSOs are highly absorbed in 
X-rays (Green et al. 1995, 1996, 2001, Gallagher et 
al. 2002). AU the BALQSOs here have N™*'' detec- 
tions or upper limits that are > 3 x lO^^cm"^. An even 
larger fraction may harbor undetected warm absorbing 
material that is too smooth in its velocity and/or geo- 
metric distribution to show distinct absorb ing troughs 
(|Gierhhski fc Dond[200l iGreen et al.l[200l . The fre- 
quency of w arm X-ray absorbers has been shown to be 
about 50% (jPorquet et al.ll2004D in low-redshift (PG) 
QSOs. Again, the bias of (color- plus broad-line-based) 
optical selection decreases their number here. The frac- 
tion of detectably-large NJj^*^ in the fully X-ray-selected 
ChaMP is larger, but nearly all such examples are found 
in objects which appear optically as narrow-line AGN 
or absorption line galaxies (XBONGs), as shown by 
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FlG. 15. — Best-fit log intrinsic absorption vs. (0.5-8 keV) counts 
for 10 random exposure-time subsamples scaled to achieve 2%, 5%, 
10%, 20%, and 50% of the full original exposure time (100%, at the 
far right) for one bright X-ray source observed by Chandra. Best- 
fit values of N^'"" are shown as blue dots, with the corresponding 
90% upper confidence limit shown as black arrows. This plot shows 
a clear trend, reflecting the skewed (one-sided positive definite and 
logarithmic) nature of the Af|p*'' parameter. 
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Fig. 16.— LEFT: Best-fit intrinsic absorption vs. (0.5-8 keV) 
counts for detected QSOs. Object type symbols are as described in 
Fig[3] Green arrows indicate upper limits, while detections (where 
the 90% lower limit exceeds log N5y*''=20) are red dots with 90% 
errorbars. The strong decrease in the upper limit envelope traces 
the increased sensitivity as a function of detected counts. Only 
a handful of objects show significant absorption above 200 counts 
where simultaneous fits to F and N^*'' are performed. RIGHT: 
Best-fit intrinsic absorption vs. redshift. The strong increase in 
the upper limit envelope refiects both the decrease in counts and 
the decreased sensitivity as a function of redshift, due to the rest- 
frame absorbed region below ~2 keV redshifting out of the Chandra 
bandpass. 

[Silverman eFall (|2005a[ ). 

6. X-RAY TO OPTICAL SEDS 

Beyond the desire to know of real trends in the 
SEDs of QSOs for the sake of understanding accretion 
physics, such trends influence other scientifically impor- 
tant research. Corrections to derive bolometric lumi- 
nosities {Liioi) from measured Lx depend strongly on 
SED trends. L\,oi in turn is key to constraining such 
fundamental parameters as accretion eff iciency and /or 
SMB H densities in the Universe (e.g., [Marconi et al.l 
2004j). or active accretion lifetim es or duty cycles (e.g., 
Hopkins et an[2005[ : [Adelberger fc Stcidel 2005). As an- 
other example, to compare the AGN spa ce densities 
found via op tical vs. X-ray selection, [Silverman et al.l 
([2005bL [20081 ) and others use the measured trend in aax 
to convert the X-ray luminosity function (LF) to an op- 
tical LF. 

Numerous studies have debated the strength and origin 
of trends in the X-ray-to-optical luminosity ratio of opti- 
cally selected quasars. This ratio is herein (and typically) 
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characterized by the spectral slope Uqx- Most studies in- 
clude large samples mixing b oth targeted and serendipi- 
tous X-ray observations fe. g..lAvni fc Tan anbaumlll98S 



Wilkes et a l. 1994; Vigna li et alj l2003l: 



20051 iSteffen et al.l 120061 : IKellv et all I20n7all . All these 



studies, using widely varying techniques and sample com- 
pilations, have found a strong correlation between Uox 
and optical luminosity, and several have also suggested a 
residual evolution of with redshift or lookback time, 
even after the luminosity correlation is accounted for. 

While large sample sizes and appropriate statistical 
techniques (see § 16. ip have been sought in these stud- 
ies, the impact of sample selection effects combined with 
large intrinsic sample dispersions in luminosity may still 
dominate the observed correlati ons. These e ffects have 
been mode l ed an d simulated bv I Yuan et al.l ()1998D and 
iTang et all (|2007D among others, who retrieve seemingly 
statistically significant but artificially-induced correla- 
tions. They highlight the need for large, well-defined 
samples to alleviate this problem, and they also note the 
importance of minimizing the effects of a strong L — z 
correlation induced by a single sample flux limit. The 
ChaMP/SDSS sample is a large step towards alleviating 
these problems. 

Contamination by unrelated physical processes should 
be eliminated wherever feasible. As mentioned above, it 
is well-known that RL QSOs tend to be X-ray bright, 
and there is evidence that a distinct, jet-related physi- 
cal process produces that extra X-ray emission. BALQ- 
SOs, on the other hand, tend to be X-ray weak, because 
of intrinsic absorption from the warm ( ionized) winds 
(jGreen et al.l 120011: [Gallagher et al.ll2002D . When inves- 
tigating distributions and correlations amongst Uoxi lo, 
and Ix, most previous authors have chosen to eliminate 
both BAL and RL QSO, to minimize contamination of 
(presumably) pure accretion-dominated X-ray emission. 
It is important to note that this precaution is neither 
complete, nor perhaps correct. For example, the intrin- 
sic radio-loudness distribution of quasars does appear to 
be well-modeled by a qua si-normal distribu tion with a 
5% tail of RL objects (Cira suolo et al.|[2003D . but the re- 
lationship of that distribution to X-ray emission is not 
well-characterized. Even after removal of quasars with a 
radio/optical flux ratio above some limit, radio-related 
X-ray emission may pervade quasar samples, and affect 
the distributions we study. Similarly, the best quanti- 
tati ve measures of absorp tion profiles via e.g., 'balnic- 
ity' (jWevmann et al.lll99lD indicate that the fraction of 
quasars with intrinsic outflows ma y be significantly un- 
derestimated (jReichard et al.|[200 3^. with classic BALQ- 
SOs just the tip of the iceberg. Equally important may 
be that BALs to date are mainly detected only beyond 
z ~ 1.4 in ground-based optical spectroscopic quasar 
samples (z ^ 1.5 in SDSS spectra), so the low redshift, 
low luminosity end of the quasar distribution may har- 
bor significant undetected warm absorption. We there- 
fore examine a variety of subsamples within and without 
these classes. When we compare bivariate regression re- 
sults I °, hkeV, and Uox, directly to Steffen et al. 
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(2006; S06 hereafter), we follow their convention of ex- 
cluding known RL QSOs and BAL QSOs, and expand 
somewhat to also eliminate QSOs with evident NALs, 
and NLSls as well. We note that their aox definition is 



the negative of the more conventional one we adopt here. 



6.1. Statistical Tools and Methods 

Many previous studies have emphasized the impor- 
tance of including X-ray upper limits in the samples, 
which require appropriate statistical treatment using 
e.g., survival analysis. Results from survival analysis 
can depend strongly on the number of detections, and 
therefore on the detection threshold, since in fiux-limited 
samples, most detections are near that threshold. The 
Chandra fields used here span a wide range of effective 
exposure times, substantially alleviating the powerful ef- 
fects that a single survey flux limit can have (e.g., by 
creating a very tight Lx — z correlation). Previous stud- 
ies have managed to limit the fraction of upper lim- 
its by using most ly targeted X-ray observations (e.g., 
iWilkes et aLlll994D . where exposure times are often cho- 
sen based on o ptical flux, or by jud icious choice of sub- 
samples (e.g., IStrateva et al] 120051: ISteffen et ahl 120061) 
desi gned to fill the lumin osity-redshift plane. For exam- 
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pie. IStrateva et all 12005) uses SDSS spectroscopic AGN 
from the DR2, observed by ROSAT PSPC for more than 
llksec. For this purpose, the DR2 AGN sample is al- 
ready biased, since many spectroscopic targets are cho- 
sen as X-ray ROSAT or FIRST radio detections. Other 
AGNs they include to fill in the L — z plane were specif- 
ically targeted for Chandra or XMM-Newton observa- 
tions for a va r iety o f reasons. A slightly larger study by 
ISteffen et all (|2006D added objects using mostly photo- 
metric AGN classifications and redshifts from COMBO- 
17, together with the Extended Chandra Deep Field 
(CDF) South (ECDF-S), but also several other smaU 
bright samples for the low L — z region. Even after se- 
lecting and combining various samples with high detec- 
tion fractions, the use of Survival Analysis techniques to 
incorporate limits must sometimes be abandoned. In 
summary, if significant selection biases may affect either 
the constituent subsamples or their ensemble, neither the 
inclusion of upper limits, nor the use of complex statisti- 
cal analysis methods should convey the impression that 
statistical results are as robust as those from a uniform, 
complete, and well-characterized sample. 

It is also worthwhile to consider the fact that some pre- 
vious studies have appended low-luminosity subsamples 
(Seyferts Is). Our SDSS sample excludes some of these 
objects because they would be spatially resolved. This 
may bias our sample in the sense that for low-redshift 
AGN (e.g., 107 objects in the MainDet sample with 
z < 0.55) we include only those that are compact (opti- 
cal light distribution consistent with the expected SDSS 
PSF). Most previous studies have instead attempted to 
include such objects, but then they correct for the inclu- 
sion of substantial host galaxy emission via e .g., spectro- 
scopic template fitting (jStrateva et al.l [20051 ). This may 
bias those samples in a different way by excluding the 
host contribution only for nearby objects. 

6.1.1. Univariate Analyses 



■^■^ For example. ISteffen et all 1)20061 ) simply drop the X-ray limits 
when regressing with I2 knV the independent variable. We 
choose instead to incorpor ate the limits u sing the two-dimensional 
Kaplan-Meier (2KM) test l lSchmitUll985l) . 
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To compare two sample distributions, we test the 'nuU' 
hypothesis that two independent random samples sub- 
ject to censoring (e.g., T for RL and RQ QSO) are ran- 
domly drawn from the same parent population. The pro- 
grams for two sample tests that we use are the Gehan, 
or generalized Wilcoxon test, and the logrank test in 
ASURV. Again, we require Pmax<0-05 to call the dis- 
tinction significant. For Pmaa;<0.10, we call the differ- 
ence "marginal" . The Kaplan-Meier method we employ 
to estimate the mean of a distribution allows for the in- 
clusion of censored values (upper or lower limits). 

6.1.2. Bivariate Regressions 

Except where noted, all the correlations studied be- 
tween X-ray and optical luminosity, and between aox 
and lo are highly significant: P < 10"^ by Cox Pro- 
portional Hazard, Kendall's r or Spearman's p tests, as 
implemented in the ASURV (Survival A nalysis for As- 
tronomy) package (jLaVallev et al.l I1992D . We deem a 
correlation significant if the maximum probability Pmax 
from all three tests is 0.05 or less. We perform bivari- 
ate linear regression s using the two -dimensional Kaplan- 
Meier (2KM) test (|Schmittl Il985f) as implemented in 
ASURV, which permits linear regression with limits in 
either axis.^^ We use 20 bootstrap iterations for error 
analysis, (more than sufhcient given the large sample size 
here) and 20 bins in each axis, with origins 27.3 (23.0, 
1.0) for I ^ ° ( hkeV, o:ox), except where samples have 
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been explicitly restricted in luminosity, or by object type 
to have N < 200, whereupon we use 10 bins. 

For bivariate regressions between X-ray and optical lu- 
minosity, we make no assumptions about which luminos- 
ity constitutes the dependent or independent variable, 
and calculate the mean (bisector) of the ordinary least 
squares (OLS) regression lines which minimize residuals 
for Y{X), X{Y), respectively.^^ While the intrinsic de- 
pendence of aox{ lo) or Uoxi Ix) is unknown, most stud- 
ies assume aox to be the dependent variable, and quote 
slopes accordingly, so we will follow suit largely for pur- 
poses of comparison. We further caution that in samples 
with large dispersions, different regression methods can 
yield different resul ts (see a rece nt mathematical review 
of related issues bv lKeIWI[2007bh . 

For completeness, we provide statistical results in table 
format for a variety of subsamples which may be of inter- 
est, even beyond those discussed in the text. For ease of 
reference, tables listing bivariate statistical results (Ta- 
bles [5] - [7]) list samples in the same order as the table 
defining samples (Table [1] 

6.2. aox for QSO Subsamples 

Here we use the D2L sample, with high detection frac- 
tion, to examine the aox distributions of several subsam- 
ples, including targets. The mean aox for 1269 QSOs 
in the D2L sample is 1.429±0.005 with median 1.370. 
Means and medians for the D2L sample and for subsam- 
ples discussed in this section are listed in Table [H as are 
result for two-sample tests. 

Our results from the other 2 bivariate regression algorithms 
in ASURV (the Buckley-James method and the parametric EM 
Algorithm) are quite consistent. 

^® In the absence of limits, these resul ts reduce reliably to the 
bisectors found by the SLOPES program IIFeieelson fc BaHulll992l : 
IBabu fc Feigelson|[T99^ : llsobe et al.lll990h . 



For 31 RLQSOs in the D2L sample, the mean 
ao^=1.377±0.028 with medi an 1. 392. Using the two- 
sample tests described in § 16.1.11 above, this distribu- 
tion is only marginally different than for the 1238 non- 
RLQSOs in the D2L sample. RL QSOs are thought to 
be more X-ray loud than RQ QSOs. 

For the 23 known BALQSOs in the D2L sample, the 
mean aox is 1.717±0.028 with median 1.66, and the dis- 
tribution is significantly different than for the non-BAL 
the D2L sample. The apparent X-ray weakness has been 
shown to be consistent with intrins ic absorption of a 
normal underlying X -ray continuum (jOreen "erani200l 
IGallaeher et"al1[200l ). 

For the 8 known NALs in the D2L sample, the mean 
ao2;is 1.463±0.056 with median 1.5, but the poor statis- 
tics render the distribution indistinguishable from the 
overall D2L sample. 

For the 19 known NLSls in the D2L sample, the mean 
ao2;is 1.54±0.08 with median 1.43, again indistinguish- 
able from the overall D2L sample. 

6.3. X-ray luminosity I2 keV vs. optical I 

2500 -/i 

Fig [17] shows a highly significant correlation of I2 keV 

with / ° , and plots our best-fit regression lines. The 
2500 A ^ ^ 

bisector regression relationship for the SDSS/ChaMP 
sample (the Main sample, limits included, 2308 QSOs), 
is 

\og[h kev) = (1.117 ± 0.017) W(?25oo;) " (7-59 ± 0.64) 

This slope is close to linear, and significantly different 
than the bisector slope /3 = 0.72 ±0.01 derived by Stef- 
fen et al. (2006; S06 hereafter) from their smaller, more 
diverse sample. The subsample philosophically closest 
to that of S06 is the NoRB sample, or its high detec- 
tion fraction version the D2LNoRB sample, since they 
exclude known RL and BAL QSOs (as well as NALs and 
NLS2s), but include targets. The D2LNoRB sample bi- 
sector slope is /3 = 1.115 ± 0.015. Our bisector slope 
results in Table [S] are closer to linear than S06 for all 
samples tested across the full luminosity range, includ- 
ing those that omit limits altogether. Fig [18] shows the 
MainDet sample (detections only) across a smaller lumi- 
nosity range, to highlight the different QSO t ypes. 

A n early linear result was also found by iHasingeTI 
(|2005() for a sample of Type 1 QSOs that spanned a sim- 
ilarly large range of luminosities as our own. That work 
used an X-ray selected sample with a high (^^95%) com- 
pleteness for spectroscopic identifications, and concluded 
that the non-linear trends seen in optically selected sam- 
ples probably result from selection effects. 

6.4. aox vs- Optical Luminosity 
Fig[Tn] (left) shows the trend of aox with I - for the 

2500 -/I 

Main sample. While the detected QSOs appear as a large 
'blob', this is due to a combination of the SDSS depth, 
the QSO optical luminosity function, and the highly effi- 
cient QSO selection in the corresponding redshift range 
1<2;<2.5 (see Fig[S]). There is a highly significant correla- 
tion: P < 10"'' by Cox Proportional Hazard, Kendall's r 
or Spearman's p for any of the samples considered (with 
or without limits). Our best-fit regression line minimiz- 
ing residuals only in aox (e.g., a Y{X) regression) for 
the Main sample is 
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Fig. 17. — X-ray (2keV) vs. optical (2500A) luminosity for the 
the Main sample. See Fig[3]for symbol types. The long-dashed ma- 
genta line is of slope unity, normalized to the sample mea ns. The 
flatt est fit is the best-fit (OLS bisector) relation from Steff eli~et al.l 
(|200Q) (their Eq. Ic), shown as a solid cyan line, spanning the lu- 
minosity range of their compilation. The best-fit OLS bisector 
regression line to the SDSS/ChaMP the Main sample (including 
non-detections) is shown in blue, spanning the full plot. Short- 
dashed lines plotted using l-cr statistical errors to the ChaMP fits 
are so close to the best-fit as to be barely discernible on the plot. 
The red line is a Y(X) regression on the same data, illustrating 
the sensitivity of fits to regression method in high dispersion data. 
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Fig. 18.— Zoom-in of the X-ray (2keV) vs. optical (2500A) 
luminosity plot for the MainDet sample. See Fig |3] for symbol 
types. Here we show a smaller luminosity range, and omit limits 
to highlight the classes indicated in the key. RL QSOs clearly 
populate the upper end of the distribution in X-ray luminosity, 
while BALQSOs are underluminous in X-rays. The fits shown here 
are performed on the MainDet sample with detections only. Line 
types are the same as in Fig |17l 



(0.061 ± 0.009) lo.g(?^^^^^) 



(0.319 ±0.258) 



This is significantly flatter than the results of S06 
(0.137 ± 0.008), consistent with the more closely linear 
relationship we find above between X-ray and optical lu- 
minosity. As can be seen from Fig [19] (bottom) and Ta- 



ble [71 keeping targets but removing RL and BAL QSOs 
(the D2LNoRB sample) yields similarly flat slopes. 

The classical Seyfert 1/QSO dividing line at Mb ^ 
—23 corresponds here to 

l''^?(^250oi) r'^-^ (°^l°S.(^^^°o^25poi) ^44.9). While 
this conventional partition is essentially arbitrary, it does 
represent a sharp discontinuity in the luminosity his- 
togram of the current sample: six times as many ob- 
jects have "QSO-like" optical luminosity as "Seyfert- 
like" . The hiLo sample, restricted to luminous QSOs 
as above, yields a slope for the aox{ ^2500!^ relation of 
0.128 ± 0.007, most similar to S06 and previous work. 
The SDSS / ChaMP sample does boast a larger number of 
low optical luminosity objects than most previous analy- 
ses, in large part because of the sensitivity of the Chan- 
dra observations, and these lower luminosity objects may 
be responsible in part for our different results across 
the full ^2500 1 range. Analysis of the (much smaller, 
with 7Vrfet=176 and 7Vtot=260) the D2LSyl sample with 
Seyfert-like luminosities, suggests that no significant cor- 
relation exists {P > 0.6 for all tests). This highlights that 
the trend of a„r (I °) may not be linear, as also found 

by S06 and iKellv et al.l (|2007af ). or may only apply for 
high luminosities. 

Note that we have investigated whether the details 
of our luminosity and aox calculations affect the re- 
sults. No signiflcantly different scientific conclusions re- 
sult from our use of best-flt F (which affects the X-ray 
iiT-correction) to calculate the X-ray luminosity. For ex- 
ample, the the MainDet sample aox{ lo) regression that 
results from instead using fixed F=1.9 and best-fit N^**" 
has slope 0.077 and intercept —0.872. 

Inasmuch as aox probes intrinsic accretion processes, it 
probably samples the balance between optical/UV black- 
body thermal emission from a geometrically thin but 
optically thick accretion disk, against X-ray emission 
from a hot, optically thin corona that upscatters the 
seed UV photons from the disk. Various physical mod- 
els can explain an a„T - I ° correlation with plausi- 
^ 2500 A , , ^ , 

ble parameters, e.g., a d isk truncation radius th at in- 
creases with luminosity (jSobolewska et al.l l2004l ). We 
also know that aox measurements can be affected by 
intrinsic absorption. This is proven by the extreme ex- 
ample of the BA LQSOs here and elsewhere ([Green et al.l 
l2001l : lGallagiier et al, 2002). Further evidence may come 
from a recent sample of AGN host galaxies of 5 clusters 
observed hy Chandraa.nd the ACS onboard HST, where 
iMartel et all ^007) found that the X -ray to optical flux 
ratio of the AGN correlates with the inclination angle of 
the host galaxies^^ If the obs erved trends are n ot domi- 
nated by selection effects (e.g.. lTang et al.l[2007( ). it seems 
likely that instrinsic absorption acts to increase the dis- 
persion of an intrinsic relationship which is dominated 
by accretion physics. 



7. EVOLUTION OF «o 

In a sample with a strong 

will naturally correlate strongly with redshift as well, due 
to the powerful redshift- luminosity trend shown in Fig[Sl 



,( I ° ) correlation, tto 



The X-ray/optical vs 
but not early-type galaxies 



inclination correlation holds for late- 
so may not apply directly to Type 1 
QSOs if they are mostly in elliptical hosts. 
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Fig. 19.— Ooi, vs. optical 2500A log luminosity for the Main sample (LEFT) and the D2LNoRB sample (RIGHT). See Fig [3 for symbol 
types. The best-fit OLS y{X) regression for each ChaMP sample is shown as a red line with errors. The best-fit relation from lSteffen et ahl 
II2006I ) is shown as a solid cyan line. 



To determine whether any rcdshift evolution of aox oc- 
curs independent of its luminosity dependence, we use 
two methods. 

First, we examine a subsample with a narrow range 
in I ° but a reasonably broad range in redshift. the 

2500^1 ° 

zBox sample (Table [7]) contains all the MainDet sample 
objects with 30.25<log( I ;)<31 and 0.5< z <2.5. 

2500 jt. 

This sample shows no significant correlation between a^x 
and redshift [Pmax — 19%). Accordingly, the nomi- 
nal best-fit regression has a slope consistent with zero 
(-0.001±0.014). 
Next, we can subtract the best-fit anxi I °) 

' , , ^ 2500 

regression to the more luminous hiLo sample 
(log(/ °) >29.8), where a simple linear fit seems 
applicable, and look for any residual aox{z) dependence 
(evolution). The expected aox{z) trend is significant 




Fig. 



RodahifL 

20.— LEFT: o 



vs. redshift for the zBox sample, re- 
stricted in log I ' to the range 30.25 — 31. No significant trend 
exists. See Fig|3]for symbol types. RIGHT: Aaox vs. redshift for 
the hiLo sample (log ( / ■)>29.8) after subtraction of its best- 



( / ° ) relation in Table [7] Although the redshift range 

2500 ./T. 



fit aox 

remains wide for this subsample, no trend is apparent 



(I ') and 

^ 2500A^ 



for the hiLo sample, induced by the Uox 

I "(z) relationships. We then subtract the best-fit 
2500 A ^ ' 

regression trend from Table [71 The residual Aaox 
shows no trend with redshift in Fig [20] (bottom) or in 
correlation tests {Pmax = 0.80). 

The large subsample sizes afforded by the 
SDSS/ChaMP QSO sample allow us to conclude, 
without resort to more elaborate statistical analyses, 
that any apparent evolutionary trend can be accounted 
for by the I ° (z) correlation in our sample, and 

2500 A ^ ' ^ ' 

that such evolution is at best very weak in the range 
0.5 < z < 2.5. 

7.1. aox 'VS. X-ray Luminosity 

A weaker trend has al so been noted in the correlation 
betw een aox and I2 keV ([Green et al.lll995l : ISteffen et al.l 
|2006[) . For the Main sample, we find 

aox ^ (0.003 ± 0.010) logih kev) + (1-384 ± 0.261) 

which (while the correlation is significant) is consistent 
with zero slope. With a higher detection fraction, the the 



D2L sample yields slope -0.0280±0.0087, so that objects 
more luminous in X-rays are also X-ray brighter (rela- 
tive to / ")■ Further removing RL and BAL QSOs 

2500 A' ° ^ 

(the D2LNoRB sample) does not change the best-fit pa- 
rameters. While the effects of different samples on the 
measured regression is significant, the aox{ hkev) re- 
lationship is particularly affected by limits, since they 
affect both axes. 

An apparent luminosity dependence of aox is gener- 
ated artificially if the intrinsic dispersi on in optical lumi- 
nosity gp exceed s that for X-rays ([Yuan et al.l [19981 : 
[Tang et an[2007f ).^^ The significance of the induced cor- 
relation is proportional to a'^/Al'^ where AIq is the op- 
tical/UV luminosity range of the sample. The magni- 
tude of the biases also depends on the luminosity func- 
tion and sample flux limits. Given these effects, the ap- 
parently strong correlations so far pu blished are all con - 
sistent with no intrinsic dependence ([Tang et al.l [2007^. 
The most convincing remedy is likely to be a volume- 
limited sample that spans a large range of both red- 

Any correlation between a dependent variable B which is de- 
rived via B oc from an independent variable A will be similarly 
affected in samples with large dispersion. 
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Fig. 21. — TOP: Best-fit X-ray spectral power-law slope F vs. 
Oox for the MainDet sample. The best-fit OLS regression relation 
for this sample is shown with a red line, and associated errors in 
dashed lines. Large open green circles show QSOs with more than 
200 (0.5-8 keV) counts and simultaneous F-A^ip*'' fits. See Fig[3]for 
other symbol types. BOTTOM: Mean F values are shown at the 
mean redshift for QSOs in 4 redshift bins of width Aa^x = 0.25 
for the MainDet sample. Errorbars show the error in the means 
for both axes. 

shift and luminosity with high detection fraction. This 
requires c areful treatment of large corn bined samples, 
similar to ISilverman et all ()2005bl . l2008h . including de- 
tailed pixel-by-pixel flux (and consequently volume) lim- 
its. While the groundwork has been laid by the ChaMP's 
xskycover analysis, such a project is well beyond the 
scope of this paper. 

7.2. aox vs. r 

Fig [21] shows best-fit F plotted against Uox for the 
MainDet sample. We detect for the first time a signifi- 
cant but shallow correlation between F and ■ Quasars 
that are relatively X-ray weak (larger aox) tend to have 
softer continuum slopes (larger F). Fig [51] shows the 
best-fit regression relation for the MainDet sample, 

F = (0.188 ± 0.106) aox + (1.676 ± 0.153). 

The measured slope of the correlation is likely some- 
what 

suppressed by the warm absorbers commonly found 
even in Type 1 AGN. In Fig[nil for the hard/weak re- 
gion bounded by > 1.6 and F < 1, we find a BAL 
(indicated by open black squares) is visible for every ob- 
ject for which a spectrum exists covering blucward of 
restframe CIV A1550. It is likely that most if not all ob- 
jects in this region are BALs, as could be determined 
e.g., from restframe UV spectroscopy. Again, these ob- 
jects are probably not intrinsically flat, but rather have 
a hard best-fit F due to undetected intrinsic absorption. 

Other samples shown in Table [8] omit targets, and 
RL or possibly absorbed QSOs all show steeper slopes. 
The steepest slope shown is for the HiCtNoTRB sample, 
which also includes only objects with counts>200, where 
N^*'' is fit independently from F. This further supports 
that absorption if anything flattens the apparent relation 
compared to the intrinsic relation. 



8. DISCUSSION 

In this study, we have presented the largest homo- 
geneous study to date of optically selected broad line 
quasars (from the SDSS) with sensitive X-ray flux limits 
(from Chandra; mode 2xl0~^''erg cm~^ s~^; 0.5-8 keV). 
Our large sample highlights the large dispersion in quasar 
properties that is unveiled whenever sensitive limits and 
wide sky areas are combined. 

We confirm and extend several well-known multiwave- 
length relations. The relationship between aox and 
2500A luminosity is confirmed for high luminosities 
(Mb < - 23, or log il^^^^;,)>29.8 (log [1^25001 ^^^^2>> 44.9) 

with slopes similar to those found previously (e.g., S06). 
Including a wider luminosity range inevitably produces 
a fiatter relation across a range of subsamples which 
test the effects of excluding Chandra PI targets, RL 
QSOs, and QSOs with BALs, or NALs as weU as NLSls. 
No significant '^oxil^^^^^) correlation exists for objects 



% detected) at lower luminosities. Another possibil- 
ity is that the relation simply flattens at low luminosi- 
ties, or has a high er order luminosity dependence (e.g., 
iKellv et al.l[2007al) . 

We find for the first time a significant, robust and 
rather steep dependence of X-ray continuum slope F on 
X-ray luminosity I2 keV in the sense that the spectrum 
hardens with increasing hkeV- A trend in the opposite 
sense has been re ported rec ently for AGN in the Chan- 
dra Deep Fields ([Saezl 120081) . but may be dominated by 
differences between Type 1 and Type 2 AGN. 

We also report a shallow but significant trend that F 
is harder for relatively X-ray bright (low aox) QSOs. We 
note that X-ray bright QSOs include the RL QSOs, and 
that RL QSOs also have predominantly fiatter (harder) 
F. We thus speculate that the overall trends of T{aox) 
and F( Z2 kev) both refiect an increase in the non-thermal 
emission fraction toward higher X-ray luminosities. Not 
all QSOs with a strong non-thermal X-ray emission com- 
ponent are necessarily radio loud - detectable radio loud- 
ness may pertain only to a fraction of these objects. Ra- 
dio bright phases may be short compared to the overall 
QSO lifetime and/or episodic. 

The black hole masses in AGN are 5-8 orders of mag- 
nitude larger than those in Galactic black hole X-ray 
binaries (GBH). Since for a given L/LEdd the disc tem- 
perature scales with mass as M"^/**, AGN disks are 
cooler than in GBH. The thermal accretion disk emis- 
sion component that dominates the soft X-ray emission 
(~2keV) in GBH corresponds to optical/UV (-2500A) 
emission in AGN. Similarly, the non-thermal emission 
(probably from Comptonized emission from the disk's 
corona) sampled as X-ray emission in AGN comes from 
harder (-20keV) X-rays in GBH. Type 1 QSOs may be 
analogous to Ga lactic black hole binarie s (GBH) in the 
high/soft state (jSobolewska et ahl |2008[ ) where a simi- 
lar trend is seen in F vs. a disc/Comptonization in- 
dex aoBH (analogous t o aox for QSOs) as we report 
here. iJester et al.l ( 20061) also compare the disk vs. non- 
thermal emission fraction of GBH and AGN, and find 
that AGN segregate by radio loudness similarly to GBH 
in regions where luminosity and/or non-thermal fraction 
are high. While the correlations we find lend significant 
support to these interpretations, the scatter is large, and 
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could be significantly reduced if extrinsic effects can be 
identified and corrected for. 

To better understand the intrinsic physics of accretion 
requires identifying and quantifying extrinsic effects such 
as absorption and non-thermal processes. Absorption 
may occur close to the SMBH gravitational radius, in 
the BLR, a molecular torus, surrounding starforming re- 
gions, the outer host galaxy or at intervening redshifts. 
The intrinsic absorption may be orientation-dependent, 
may evolve with redshift, and may be a function of lumi- 
nosity as well. Contributions from non-thermal processes 
certainly play a role, whether or not it is reflected in de- 
tectable radio emission, and that role may change with 
SMBH mass, spin, and environment, so consequently 
with lookback time as well. 

Let's face it - quasars are complicated, but on average 
they don't change much. If quasar SED changes with 
luminosity were large compared to the dispersion in the 
population and relatively immune to selection effects, we 
would have been using them for standard candles a long 
time ago. Large samples, uniformly observed and ana- 
lyzed, offer the greatest hope to disentangle the inter- 
twining mysteries. 
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APPENDIX 

The ChaMP has developed and implemented an xskycover pipeline which creates sensitivity maps for all ChaMP 
sky regions imaged by ACIS. This allows (1) identification of imaged-but-undetected objects (2) counts hmits for 50% 
and 90% detection completeness (3) flux sensitivity vs. sky coverage for any subset of obsids, necessary for logN-logS 
and luminosity function calculations and (4) flux upper- limits at any sky position. The basic recipe is as follows. We 
use the wavdetect detection algorithm in CIAO (Free man et al.l[2002l ) to generate threshold maps at each wavdetect 
kernel scale actually run (1, 2, 4, 8, 16, and 32 pixels). The threshold maps, computed from the local background 
intensity, determine the magnitude of the source counts necessary for a detection at each pixel with a detection 
threshold of F = 10~^ (corresponding to 1 false source per 10^ pixels). Thus, when a source is not detected at a 
given location, the threshold map value serves as an upper limit to the source counts. Nominally, a source with true 
intensity equal to this counts limit would be detected in approximately half the instances that the source is observed 
under the same conditions. To retain fidelity yet create a reasonably-sized and easily-sampled sensitivity table covering 
the full ChaMP, we first average the threshold map values in sky pixels, each 10 x 10 arcsec, whose boundaries are 
chosen to match a regular commensurate grid across the sky. The final value of this counts limit in any given sky 
pixel is interpolated from the two threshold maps computed at wavelet scales that bracket the size of the PSF at that 
location. 

Note that the identification of the threshold map value one-to-one with the counts limit is valid only for a specific 
shape of the PSF (see Equation 6 of Freeman et al. 2002). In particular, the strength of the putative source can vary 
widely based on the correspondence between the PSF size and the wavelet scale, as well as the shape of the PSF. For 
non-Gaussian PSF shapes (such as are found with Chandra at larger off-axis angles), the threshold values must be 
corrected before a source counts limit is determined. We calibrate this correction factor by comparing the detection 
threshold map values with simu lated source retrie val experiments conducted on a subsample of ChaMP fields - the 
130 Cycle 1-2 obsids studied by Ki m et al] (|2007bf ). While the threshold values give us a reliable map of variation on 
the sky, we must find a normalization from these simulations. Summing over a large number (^50,000) of simulated 

The 39% radius is determined using the PSF enclosed counts fraction, which corresponds to the la 2-dimensional PSF size, from the 
hbrary generated by the CXC Calibration group as documented at |htt p : / / cxc . harvard .edu/cal/Hrma/psf 
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sources, we compare the actual detection fraction from the simulations to the ratio of input (simulated) counts to 
derived (threshold) counts. From the fine binning available in these data, we interpolate to find the normalization 
that yields the correct counts values for 50% and 90% completeness, across a range of exposures, background levels, 
chip types, and off-axis angles. Our experiments show that the only significant dependence of the correction factor 
(normalization) is on off-axis angle, and that dependence warrants only a linear correction dependence with best fit 
A^50 = 1-32 + 0.198xOAA {Ngo = 2.08 -I- 0.33 1 xOAA) for 50 % (90%) source detection probability. 

Our method has been verified recently bv I Alder oft et al.| (|2008) using the Chandra Deep Field South (CDF-S) 
data available from the CXC. These data include the full 1.8Msec from 2000 (jGiacconi et al.|[2002[ ). and Director's 
Discretionary Time observations in 2007. The full 1.8 Msec exposure was used to generate a source list which extended 
to faint fluxes. The 21 individual obsids comprising the full exposure were then treated as realizations of an experiment 
to detect these sources and the ensemble of detection statistics were analyzed by the same method as used for the 
^50,000 simulated sources. We found excellent agreement for the 50% detection threshold and a slight disparity ( 10%) 
for the 90% detection threshold which could be explained by source variability. 

REFERENCES 



Adelberger, K. L., Steidel, C. C, ApJ, 629, 636 
Aldcroft, T.L. 2006, Bulletin of the American Astronomical 
Society, 38, 376 

Aldcroft, T.L., Green, P.J., Connelly, J., & Kashyap, V., & Kim, 
D.-W. 2008, Bulletin of the American Astronomical Society, 40, 
in press 

Antonucci, R. 1993, ARA&A, 31, 473 
Avni, Y., & Tananbaum, H. 1982, ApJ, 262, L17 
Babu, G. J. and Feigelson, E. D., 1992, Simulation & Computation, 
21, 533 

Ballantyne, D. R., Everett, J. E., & Murray, N. 2006, ApJ, 639, 
740 

Barger, A. J., Cowie, L. L., Mushotzky, R. F., Yang, Y., Wang, 
W.-H., Steffen, A. T., & Capak, P. 2005, AJ, 129, 578 

Barvainis, R., Lehar, J., Birkinshaw, M., Falcke, H., & Blundell, 
K. M. 2005, ApJ, 618, 108 

Becker, R. H., White, R. L., & Helfand, D. J. 1995, ApJ, 450, 559 

Blackburne, J. A., Pooley, D., & Rappaport, S. 2006, ApJ, 640, 
569 

Blundell, K. M., & Beasley A. J. 1998, MNRAS, 299, 165 
Boiler, T., Brandt, W. N., & Fink, H. 1996, A&A, 305, 53 
Boroson, T. A., & Green, R. F. 1992, ApJS, 80, 109 
Boroson, T. A. 2002, ApJ, 565, 78 

Boyle, B. J., Shanks, T., & Peterson, B. A. 1988, MNRAS, 235, 
935 

Brandt, W. N., & Hasinger, G. 2005, ARA&A, 43, 827 
Cirasuolo, M., Celotti, A., Magliocchetti, M., & Danese, L. 2003, 

MNRAS, 346, 447 
Condon, J. J., Cotton, W. D., Greisen, E. W., Yin, Q. F., Perley, 

R. A., Taylor, G. B., & Broderick, J. J. 1998, AJ, 115, 1693 
Croom, S. M., Smith, R. J., Boyle, B. J., Shanks, T., Miller, L., 

Outram, P. J., & Loaring, N. S. 2004, MNRAS, 349, 1397 
Cash, W. 1979, ApJ, 228, 939 

Czerny, B., Siemiginowska, A., Janiuk, A., & Gupta, A. C. 2008, 

MNRAS, 386, 1557 
Dai, X., Chartas, G., Eracleous, M., & Garmire, G. P. 2004, ApJ, 

605, 45 

Dai, X., Shankar, F., & Sivakoff, G. R. 2008, ApJ, 672, 108 

Dickey J.M., & Lockman, F.J. 1990, ARA&A, 28, 215 

Fabian, A. C, Iwasawa, K., Reynolds, C. S., & Young, A. J. 2000, 

PASP, 112, 1145 
Fabbiano, G., et al. 2007, Astronomical Data Analysis Software 

and Systems XVI, 376, 172 
Falcke, H., Kording, E., & Markoff, S. 2004, A&Ap, 414, 895 
Feigelson, E. D. & Babu, G. J., 1992, ApJ, 397, 55 
Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 
Freeman, P. E., Kashyap, V., Rosner, R., & Lamb, D. Q. 2002, 

ApJS, 138, 185 
Gebhardt, K., et al. 2000, ApJ, 539, L13 
Galbiati, E., et al. 2005, A&A, 430, 927 
Gallagher, S. C. et al. 2002, ApJ, 567, 37 
Gierlihski, M., & Done, C. 2004, MNRAS, 349, L7 
Giacconi, R., et al. 2002, ApJS, 139, 369 
Gilli, R., Comastri, A., & Hasinger, G. 2007, A&A, 463, 79 
Granato, G. L., et al. 2004, ApJ, 600, 580 

Green, A. R., McHardy I. M., & Lehto, H. J. 1993, MNRAS, 265, 
664 

Green, P. J., et al. 1995, ApJ, 450, 51 
Green, P. J. 1996, ApJ, 467, 61 



Green, P. J. et al 2001, ApJ, 558, 109 
Green, P. J., et al. 2002, ApJ, 571, 721 
Green, P. J., et al. 2004, ApJS, 150, 43 

Green, P. J., Infante, L., Lopez, S., Aldcroft, T. L., & Winn, J. N. 

2005, ApJ, 630, 142 
Green, P. J. 2006, ApJ, 644, 733 
Grindlay J. E., et al. 2005, ApJ, 635, 920 

Grupe, D., Wills, B. J., Leighly K. M., & Meusinger, H. 2004, AJ, 
127, 156 

Grupe, D. & Mathur, S. 2004, ApJL, 606, 41 

Hewett, P. C, Foltz, C. B., & Chaffee, F. H. 1995, AJ, 109, 1498 

Hewett, P. C, & Foltz, G. B., 2003, AJ, 125, 1784 

Hasinger G., 2005, in Merloni A., Nayakshin S., Sunyaev R., eds. 

Growing Black Holes: Accretion in a Cosmological Context. 

Springer- Verlag, Berlin, p. 418 
Hennawi, J. F., et al. 2006, AJ, 131, 1 

Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., Robertson, 

B., Volker, S. 2005, ApJ, 630, 705 
Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., Robertson, 

B., & Springel, V. 2006, ApJS, 163, 1 
Hopkins, P. F., et al. 2006, ApJ, 652, 864 
Hopkins, P. F., et al.2008, ApJS, 175, 356 

Isobe, T., Feigelson, E. D., Akritas, M. & Babu, G. J., 1990, ApJ, 
364, 104 

Ivezic, Z., et al. 2002, AJ, 124, 2364 

Jester, S., Kording, E., & Fender, R. 2006, VI Microquasar 

Workshop: Microquasars and Beyond 
Jiang, L., Fan, X., Ivezic, Z., Richards, G. T., Schneider, D. P., 

Strauss, M. A., & Kelly B. C. 2007, ApJ, 656, 680 
Just, D. W., Brandt, W. N., Shemmer, O., Steffen, A. T., 

Schneider, D. P., Chartas, G., & Garmire, G. P. 2007, ApJ, 665, 

1004 

Kelly, B. C, Bechtold, J., Siemiginowska, A., Aldcroft, T., & 

Sobolewska, M. 2007, ApJ, 657, 116 
Kelly B. C. 2007, ApJ, 665, 1489 

Kelly, B. C, Bechtold, J., Trump, J. R., Vestergaard, M., & 

Siemiginowska, A. 2008, ApJ, in press, (|arXiv:0801.23839 
Kim, D.-W., et al. 2004, ApJS, 150, 19 
Kim, M., et al. 2007a, ApJS, 169, 401 

Kim, M., Wilkes, B. J., Kim, D.-W., Green, P. J., Barkhouse, 
W. A., Lee, M. G., Silverman, J. D., & Tananbaum, H. D. 2007b, 
ApJ, 659, 29 

Knigge, C, Scaringi, S., Goad, M. R., & Cottis, C. E. 2008, 

MNRASin press, ( arXiv:0802.3697 ) 
La Franca, F., et al. 2005, ApJ, 635, 864 

LaValley M., Isobe, T., & Feigelson, E. D. 1992, ADASS, 1, 245 
Lamer, G., Schwope, A., Wisotzki, L., & Christensen, L. 2006, 

A&A, 454, 493 
Lehmer, B. D., et al. 2007, ApJ, 657, 681 

Maccarone, T. J., Gallo, E., & Fender, R. 2003, MNRAS, 345, L19 
Mainieri, V., et al. 2007, ApJS, 172, 368 

Marconi, A., Risaliti, G., Gilli, R., Hunt, L. K., Maiolino, R., & 

Salvati, M. 2004, MNRAS, 351, 16 
Martel, A. R., Menanteau, F., Tozzi, P., Ford, H. C, & Infante, L. 

2007, ApJS, 168, 19 
Misawa, T., Eracleous, M., Chartas, G., & Charlton, J. C. 2008, 

ApJ, 677, 863 



Properties of SDSS Quasars in the ChaMP 



21 



Myers, A. D., Outram, P. J., Shanks, T., Boyle, B. J., Groom, 
S. M., Loaring, N. S., Miller, L., & Smith, R. J. 2003, MNRAS, 
342, 467 

Myers, A. D., et al. 2006, ApJ, 638, 622 

Myers, A. D., et al. 2007, ApJ, 658, 99 

Myers, A. D., et al. 2008, ApJ, in press, ( arXiv:0709.3474 ) 
Middleton, M., Done, C, & GierUhski, M. 2007, MNRAS, 381, 
1426 

Oke, J. B. & Gunn, J. E. 1983, ApJ, 266, 713 
Osterbrock, D. E., & Pogge, R. W. 1985, ApJ, 297, 166 
Page, M. J., et al. 2003, AN, 324, 101 

Page, K. L., Reeves, J. N., O'Brien, P. T., Turner, M. J. L., & 

Worrall, D. M. 2004, MNRAS, 353, 133 
Page, G. G. 2006, Astronomical Data Analysis Software and 

Systems XV, 351, 85 
Perola, G. G., et al. 2004, A&A, 421, 491 

Porquet, D., Reeves, J. N., O'Brien, P., & Brinkmann, W. 2004, 
A&A, 422, 85 

Ptak, A., & Griffiths, R. 2003, Astronomical Data Analysis 

Software and Systems XII, 295, 465 
Reeves, J. N., & Turner, M. J. L. 2000, MNRAS, 316, 234 
Reichard, T. A., et al. 2003, AJ, 126, 2594 
Reynolds, G. S., & Nowak, M. A. 2003, Phys. Rep., 377, 389 
Richards, G. T., et al. 2001, AJ, 121, 2308 
Richards, G. T., et al. 2002, AJ, 123, 2945 
Richards, G. T., et al. 2004, ApJS, 155, 257 
Richards, G. T., et al. 2005, MNRAS, 360, 839 
Richards, G. T., et al. 2006a, AJ, 131, 2766 
Richards, G. T., et al. 2006b, ApJS, 166, 470 
Richards, G. T., et al. 2008, ApJS, submitted 
Saez, C. et al. 2008, AJ, 135, 1505 

Scannapieco, E., Silk, J., & Bouwens, R. 2005, ApJ, 635, L13 
Schmitt, J. H. 1985, ApJ, 293, 178 

Schneider, D. P., Schmidt, M., & Gunn, J. E. 1994, AJ, 107, 1245 
Schneider, D. P., et al. 2007, AJ, 134, 102 

Shastri, P., Wilkes, B. J., Elvis, M., & McDowell, J. 1993, ApJ, 
410, 29 

Shemmer, O., Brandt, W. N., Vignali, C., Schneider, D. P., Fan, 

X., Richards, G. T., & Strauss, M. A. 2005, ApJ, 630, 729 
Shemmer, O., Brandt, W. N., Netzer, H., Maiolino, R., & Kaspi, 

S. 2006, ApJ, 646, L29 
Shemmer, O., Brandt, W. N., Netzer, H., Maiolino, R., & Kaspi, 

S. 2008, (arXiv:0804.0803) 
Shen, S., White, S. D. M., Mo, H. J., Voges, W., Kauffmann, G., 

Tremonti, C., & Anderson, S. F. 2006, MNRAS, 369, 1639 
Siemiginowska, A., Cheung, G. C., LaMassa, S., Burke, D. J., 

Aldcroft, T. L., Bechtold, J., Elvis, M., & Worrall, D. M. 2005, 

ApJ, 632, 110 



Silverman, J. D., et al. 2005a ApJ, 618, 123 
Silverman, J. D., et al. 2005b, ApJ, 624, 630 
Silverman, J. D., et al. 2008, ApJ, 679, 118 

Sobolewska, M. A., Siemiginowska, A., &Zdotycki, P. T. 2004, ApJ, 
617, 102 

Steffen, A. T., Strateva, I., Brandt, W. N., Alexander, D. M., 
Koekemoer, A. M., Lehmer, B. D., Schneider, D. P., &; Vignali, 
C. 2006, AJ, 131, 2826 

Sobolewska, M. A., Gierlihski, M., & Siemiginowska, A. MNRAS, 
submitted Strateva, I., Brandt, W. N., Alexander, D. M., 
Koekemoer, A. M., Lehmer, B. D., Schneider, D. P., & Vignali, 

C. 2006, AJ, 131, 2826 

Steffen, A. T., et al. 2007, ApJL, 667, L25 

Strateva, I. V., Brandt, W. N., Schneider, D. P., Vanden Berk, 

D. G., & Vignali, C. 2005, AJ, 130, 387 
Tajer, M., et al. 2007, A&A, 467, 73 

Tang, S., Zhang, S. N., & Hopkins, P. F. 2007, MNRAS, 377, 1113 
Treister, E., & Urry, C. M. 2006, ApJ, 652, L79 
Trump, J. R., et al. 2006, ApJS, 165, 1 

Ueda, Y., Akiyama, M., Ohta, K., & Takamitsu, M. 2003, ApJ, 
598, 886 

Vignali, C., Brandt, W. N., & Schneider, D. P. 2003, AJ, 125, 433 
Vignali, C., Brandt, W. N., Schneider, D. P., & Kaspi, S. 2005, AJ, 
129, 2519 

Warner, C., Hamann, F., & Dietrich, M. 2004, ApJ, 608, 136 
Weinstein, M. A., et al. 2004, ApJS, 155, 243 

Weymann, R. J., Morris, S. L., Foltz, G. B., & Hewett, P. G. 1991, 
ApJ, 373, 23 

White, R. L., Helfand, D. J., Becker, R. H., Ghkman, E., & de 

Vries, W. 2007, ApJ, 654, 99 
Wilkes, B. J., Tananbaum, H., Worrall, D. M., Avni, Y., Gey, M. S., 

& Flanagan, J. 1994, ApJS, 92, 53 
Williams, R. J., Mathur, S., & Pogge, R. W. 2004, ApJ, 610, 737 
Wilman, R. J. & Fabian, A. C. 1999, ApJ, 522, 157 
Wolf, G., Wisotzki, L., Borch, A., Dye, S., Kleinheinrich, M., & 

Meisenheimer, K. 2003, A&A, 408, 499 
Worrall, D. M., Tananbaum, H., Giommi, P., & Zamorani, G. 1987, 

ApJ, 313, 596 
Wyithe, J. S. B., & Loeb, A. 2002, ApJ, 581, 886 
Wyithe, J. S. B., & Loeb, A. 2005, ApJ, 634, 910 
York, D. G., et al. 2000, AJ, 120, 1579 
Yu, Q., & Tremaine, S. 2002, MNRAS, 335, 965 
Yuan, W., Brinkmann, W., Siebert, J., & Voges, W. 1998, A&A, 

330, 108 

Zamorani, G., et al. 1981, ApJ, 245, 357 



TABLE 1 
Quasar Sample Definitions 



Sample 


Lims 


Targets 


RL 




OAA Naet 




^total 


%Det 










Primary Samples 










Main 


y 








1135 


1173 


2308 


49 


MainDet 










1135 





1135 


100 


noTDet 




n 






1053 





1053 


100 


D2L 


y 






4 


12 922 


347 


1269 


72 


D2LNoRB 


y 




n 


n 4 


12 866 


338 


1204 


71 


hiLo 


y 




n 


n 


847 


961 


1808 


46 


HiCtNoTRB 




n 


n 




129 





129 


100 










Other Samples 










NoRB 


y 




n 


n 


1054 


1144 


2198 


47 


NoRBDet 






n 


n 


1054 





1054 


100 


D2LNoTRB 


y 


11 


n 


n 4 


12 828 


338 


1166 


71 


hiLoLx 






n 


n 


801 





801 


100 


zLxBox 










817 





817 


100 


LoBox 










530 





530 


100 


zBox 


y 








420 


360 


780 


53 


zBoxDet 










420 





420 


100 


D2LSyl 


y 




n 


n 


176 


84 


260 


68 


HiCt 










157 





157 


100 


HiCtNoTRB 




n 


n 


n 


129 





129 


100 



" If 'y', sample includes X-ray non-detections. 

^ If 'n', sample excludes QSOs with evident BALs or NALs and also NLSls. 



Table 2. PROPERTIES OF SDSS QUASARS DETECTED BY CHANDRA 
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,733 





.145 


0, 


,933 





.140 





.240 





.2053 


s 


J002752. 4+002615 


XS04080B7. 


-003 


4.5 


39.2 


7, 


,6 


1 


.5 


44.9 


5880 15 509 809G 599 37 


7. 


,88101 


0, 


,572282 


18, 


,474 


1 


.875 


0, 


,868 


1 


.620 


2 


.040 


1 


.7354 


s 


.7003131.4+003420 


XS02101B7. 


.002 


1.2 


28.8 


6. 




3 


.8 


2.4 


587727227305066749 


10. 


.07393 


-9, 


,190477 


20 


,752 





.145 


0, 


,649 





.100 





.250 





.1450 




.7004017.7-091125 


XS04888B3. 


.010 


7.6 


19.8 




,8 


8 


.8 


3.4 


587727227305197873 


10, 


.30287 


-9. 


,238581 


20, 


,770 


1, 


.295 


0, 


,635 


1. 


.000 


1 


.570 


1. 


.2950 




J004112. 6-091417 


XS04888B1. 


.010 


11.0 


16.1 


6, 


,5 


8. 


.6 


3.4 


587731185135648990 


12, 


.47624 


-0. 


,939257 


20, 


.389 


2, 


.145 


0, 


.702 


1, 


.440 


2 


.260 


2, 


.1450 




J004954.3-005620 


XS04825B7. 


.018 


4.7 


10.8 


4, 


.4 


12, 


.9 


2.1 


587731185135648996 


12, 


.48751 


-0. 


,968442 


20, 


,758 





.385 


0, 


.551 





.250 





.500 





.3850 




J004957.0-005806 


XS04825B7. 


-012 


5.8 


26.9 


6, 


.6 


12 


.9 


2.1 


588015508201144501 


12. 


,62757 


-0, 


,780014 


20 


,845 


1 


.265 


0, 


,961 





.900 


1 


.470 


1 


.2650 




.70050.30.6-004649 


XS04825B2. 


-003 


8.5 


21.9 


6, 


,2 


12 


.3 


2.1 


588015508201144513 


12. 


,65482 


-0, 


.808028 


20, 


,411 


1 


.855 


0, 


,883 


1 


.620 


2 


.110 


1 


.8550 




.7005037.2-004829 


XS04825B2. 


.001 


8.8 


23.8 


6, 


,2 


12 


.3 


2.1 


587731186209783863 


13. 


.47981 


-0. 


,052600 


17, 


,984 





.485 


0, 


,778 


0. 


.390 





.700 


1. 


.7189 


s 


J005355. 1-000309 


XS04830B7. 


.001 


5^ 


27.7 


6, 


,4 


7. 


.0 


1.9 


588015509275803698 


14, 


.77296 


0. 


,114358 


17, 


.489 


0, 


.745 


0, 


,842 


0, 


.660 





.950 


0, 


.7189 


s 


J005905.4+000651 


XS02179B6. 


-001 


324.7 


19, 


.6 


2, 


.2 


3.0 


588015509275869378 


14, 


.84438 


0. 


,050395 


19, 


.194 


4, 


.385 


0, 


.998 


4, 


.200 


4 


.560 


4, 


.1544 


s 


J005922.6+000301 


XS02179B7. 


-003 


to 


10.8 


4, 


.4 


2, 


.6 


3.0 



Table 2 — Continued 



SDSS ObjID N'ff*^ N^' JV|^ r r*** r'° log/^ log hkeV log '2500 A """^ f2acm Ext fi Class Targ Comments 



(1) 


(18) 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 




(25) 


(26) 


(27) 


(28) 


(29) 


(20) 


(31) 


(32) 


(33) 


ocS77oiio7z77ooyDyo 








1, 


.53 


0, 


,45 


-0, 


,44 


-12, 


.960 


26 


,988 


30, 


,590 


1 


,383 







1 


,06 










COTTOl 1 TTri CI C AO 

oof 1 d 1 Iti i Z 1 lybbUoo 








2 


,40 


0, 


,42 


-0, 


,38 


-12, 


,850 


27 


,065 


30 


,978 


1 


,502 










,62 





1 




5880 155 10343385 196 








2 


,14 


1 


,25 


-1 


,03 


-13, 


,990 


26 


,286 


30 


,391 


1 


,576 







1 


,54 










oof f 31 18DzU4D0ObDD 


0.00 


0.27 


0.00 


1 


,74 


0, 


,30 


-0, 


,24 


-12, 


321 


27 


,308 


30 


,020 


1 


,041 


3oy ( .6U 


1 


4 


,86 








iNhjU: h hJ^Jb .J00UO-UUU4 


00/ f 31 18dzU4dOO f 04 








1 


,31 


0, 


,49 


-0, 


,46 


-12, 


,870 


27 


,228 


29 


,899 


1 


,025 







1 


,75 










OooUlOoUo 1 33^311 1 1 








2, 


,05 


0, 


,46 


-0, 


,44 


-13, 


,180 


26 


,992 


31, 


,173 


1 


,605 







0, 


,59 





1 


iNiLU: Lihil^D 0UU4-UU3J 


588015508733231262 








3, 


,50 


1 


,87 


-1, 


,26 


-14, 


,000 


26 


.563 


30, 


,650 


1 


,669 







1 


,68 










588015508733231265 








2, 


,27 


0, 


,68 


-0, 


,62 


-13, 


,470 


26 


.850 


30, 


,758 


1 


,500 







1 


,21 










587731186204737774 








1 


82 


0, 


,91 


-0 


,79 


-13, 


230 


26 


.561 


29 


,765 


1 


,230 







1 


,73 










5S7730773889974538 








1 


88 





,37 


-0 


35 


-12 


580 


28 


.701 


32 


,063 


1 


,290 










,75 





1 




587731 186742198291 








2 


,13 


1 


,02 


-0 


89 


-13 


730 


26 


.681 


31 


,020 


1 


,666 







1 


,03 










588015509270822924 








2 


,10 


0, 


,87 


-0 


74 


-13 


820 


26 


.587 


31 


,220 


1 


,778 










,82 


1 


1 


HiBAL NED: LBQS 0010-0012 


588015509270823186 








1, 


,47 


0, 


,43 


-0, 


,41 


-13, 


220 


27 


. 142 


30, 


,316 


1 


,218 







1 


,69 










588290881639350481 








2, 


,05 


1, 


,48 


-1, 


,24 


-14, 


,640 


24 


.819 


29, 


,639 


1 


,850 







1 


,69 










588290881639350569 








1. 


,72 





,73 


-0. 


,67 


-14, 


,010 


26 


.412 


30 


,241 


1 


,470 







1 


,82 










587730775501504810 








1 


69 





,70 


-0 


63 


-13, 


,450 


24 


.354 


27 


,578 


1 


,238 







1 


,87 










588290881639350397 


0.00 


0.47 


0.00 


1 


77 





,16 


-0 


12 


-12, 


,609 


27 


.729 


31 


,676 


1 


,515 










,31 





1 




588015507661324390 








1 


,94 


0, 


,60 


-0, 


,57 


-13, 


,000 


26, 


.799 


30, 


,231 


1 


,317 







1 


,31 










588015509809266720 








1, 


,28 


0, 


,50 


-0, 


,46 


-12, 


,540 


25, 


.418 


29, 


,033 


1 


,388 







0, 


,53 


3 


1 


NLSyl 


588015509809659937 








2, 


,23 


0, 


,58 


-0, 


,54 


-13, 


,360 


26, 


.822 


30, 


,939 


1 


,681 







0, 


,83 










587727227305066749 








1, 


,91 


0, 


,73 


-0, 


,67 


-13, 


,540 


24, 


.084 


27, 


.718 


1 


.395 







1 


,74 










587727227305197873 








1, 


,16 


0, 


,79 


-0, 


,73 


-13, 


,410 


26, 


.458 


29, 


.783 


1, 


.276 







1 


,75 










587731185135648990 




dd 
H 




0, 


,95 


0, 


,89 


-0, 


,86 


-13. 


,950 


26 


.457 


30 


.449 


1 


.532 







1 


,60 










587731185135648996 






2 


,30 





,65 


-0 


,59 


-13 


810 


24 


.775 


28 


.568 


1 


.456 







1 


,74 










588015508201144501 






2, 


,79 


0, 


,90 


-0, 


,80 


-13. 


,690 


26, 


.153 


29, 


.830 


1 


.412 







1 


,78 










588015508201144513 






1, 


,77 


0, 


,71 


-0, 


,62 


-13. 


,580 


26 


.673 


30 


.307 


1 


.395 







1 


.60 










587731186209783863 




w 

0.47 




0, 


,90 


0, 


,48 


-0, 


,48 


-13. 


,300 


26, 


.871 


31, 


.118 


1 


.630 







0, 


.63 


2 


1 


many strong NALs 


588015509275803698 


0.00 


0.00 


1, 


,61 


0, 


,26 


-0, 


,16 


-11. 


,751 


27, 


.486 


30, 


.690 


1 


.230 


2508.80 


1 


3, 


.84 










588015509275869378 








1, 


,32 


0, 


,84 


-0, 


,77 


-13. 


,410 


27, 


.688 


31, 


.642 


1 


.480 







1 


.12 





1 





Note. — The complete version of this tabic is in the electronic edition of the Journal. The printed edition contains only a sample. (1) SDSS Object ID (2) SDSS R.A. (J2000) (3) SDSS Dec (J2000) (4) SDSS 
asinh mag_psf i, dercddened. (5) Photometric rcdshift (see Weinstcin et al. 2004). (6) Photometric rcdshift range probability. (7) Lower limit of photometric redshift range. (8) Upper limit of photometric 
redshift range. (9) Best redshift: spectroscopic if different than ^p/^of . (10) Reference for spectroscopic redshift: S: SDSS, O: ChaMP, R: published reference from NED. (11) ChaMP lATJ source name. (12) 
ChaMP internal source ID, format XSoooooBc_nnn where ooooo is Chandra obsid, c is ACIS CCD ID, and nnn is source ID on that CCD. (13) Chandra Off-axis angle in arcmin. (14) net 0.3-8 keV source 
counts. (15) rms uncertainty on net counts. (16) Vignetting-corrected exposure time in ksec. (17) Galactic column in units lO'^'^ cm — 2 (18) Best-fit X-ray intrinsic column in lO'^'^ cm — 2, only included four 
counts>200. (19) 90% upper limit on intrinsic column in 10^^ cm-2 (20) 90% lower limit on intrinsic column in 10^^ cm-2 (21) Best-fit X-ray power-law index T (22) 90% upper limit on T (23) 90% lower 
limit on T (24) log X-ray flux (0.5-8 keV) in erg cm~^ s~^ (25) log X-ray luminosity at 2 keV in erg s~^ Hz~^ (26) log optical/UV luminosity at 2500A in erg s""*^ Hz~^ (27) olqx, the optical/UV to X-ray 
spectral index. (28) 20cm radio flux in mjy from FIRST or NVSS. (29) radio extent flag (30) radio loudness R (31) Spectral Class: 1= BAL, 2= NAL, 3= NLSl. (32) 1= intended Chandra PI target. (33) 
Comments 



Table 3. PROPERTIES OF SDSS QUASARS WITH CHANDRA LIMITS 



SDSS ObjID 


R.A. (J2000) Dec 




^phot 








best 


Spec Ref 


obsid 


ccdid 


OAA c 


:ounts < 


Exposure 


H 




2 kcV 


loe / 

^ 2500 A 






20cm 


R 


Class 


(1) 


(2) 


(3) 








( ') 


(S) 


(9) 


(10) 


(11) 


(12) 


(13} 


(14) 


(15) 


(16) 


(17) 


(18) 


(10) 


(20) 


(21) 


(22) 


(23) 


5 8 8 1 .'i 5 1 034338 5 2 9 5 


0.56337 


0.009429 


20.056 


4-6050 


0.949 


4.440 


5-22 


4.6050 




4861 


7 


6.4 


16.0 


4.5 


2.5 


-13.650 


27-554 


30.856 




268 




1.82 





588015510343385302 


0.63485 


0.870704 


20.819 


2-2250 


0.415 


1.840 


2-70 


2.2250 




4861 


7 


2.8 


10.0 


4.6 


2.5 


-13.865 


26-581 


30.311 




432 




1-77 





588015509270823376 


3.32406 


0.057045 


20.084 


0-0650 


0.616 


0.060 


0-12 


0.0650 




4829 


7 


3.5 


11.1 


6.0 


5.6 


-13.926 


22-055 


26.290 




280 




1-83 





588015509270888717 


3.40998 


0. 141416 


20.457 


1.4450 


0.615 


0.960 


1.61 


1.4450 




4829 


3 


9.0 


20.2 


4.4 


5.6 


-13.239 


26-747 


30.023 




258 




1.62 





588015509270888725 


3.45158 


0.098009 


20.804 


1.4950 


0.774 


0.950 


2.22 


1.4950 




4829 


3 


10.8 


22.7 


4.6 


5.6 


-13.211 


26.811 


29.915 




191 




1.76 





587730775501504550 


5.00208 


15.852781 


19.190 


1.6850 


0.915 


1.440 


2.11 


1.7510 


s 


1595 


5 


10.7 


25.3 


5.3 


4.2 


-13.405 


26.786 


30.742 




519 




1.12 





587727180601295054 


10.04659 


-9.111831 


20.220 


2.0950 


0.653 


1.290 


2.26 


2.0950 




4888 


3 


7.5 


17.2 


7.8 


3.4 


-13.668 


26.714 


30.495 




.452 




1.53 





587727180601557048 


10.60443 


-8.983756 


19.739 


1.8750 


0.916 


1.590 


2.13 


1.8750 




4886 


2 


9.9 


22.6 


6.9 


3.6 


-13.498 


26.766 


30.586 








1.34 





587727227305394556 


10.73418 


-9.141124 


20.920 


0.1450 


0.370 


0.140 


0.25 


0.1450 




4886 


1 


2.3 


9.4 


7.9 


3.6 


-13.930 


23.694 


27.029 




280 




1.81 





587727180601622916 


10.82587 


-9.032790 


20.550 


2.8550 


0.866 


2.520 


3.26 


2.8550 




4886 





8.1 


17.7 


7.0 


3.6 


-13.605 


27.104 


30.574 




332 




1.66 





587731185135649019 


12.55619 


-0.978559 


20.734 


0-3950 


0.575 


0.240 


0.08 


0.3950 




4825 


7 


6.2 


15.1 


11.4 


2.1 


-14.075 


24.536 


28-869 




663 




1.73 





588015508201111424 


12.68671 


-0.838214 


20.077 


0-5850 


0.525 


0.420 


0.68 


0.5850 




4825 


3 


10.0 


22.7 


5.9 


2.1 


-13.317 


25.702 


20,374 




410 




1-47 





588015509275279576 


13.48778 


0.041851 


21.014 


2-3550 


0.485 


2.080 


2.73 


2.3550 




4830 


6 


5.1 


12.8 


5.0 


1.9 


-13.635 


26.871 


30.327 




327 




1.85 





587731186200849918 


13.O0874 


-0. 132789 


20.138 


4-6750 


0.958 


4.480 


4.87 


4.6750 




4830 


7 


5.7 


14.0 


6.2 


1-9 


-13.843 


27-376 


31-176 




459 




1.49 





587731511532454189 


19.71203 


-0.063115 


20.433 


2.3950 


0.849 


2.190 


2.02 


2.3950 




4963 


7 


4.1 


15.0 


36.9 


14.1 


-14.575 


25-040 


30.534 




760 




1.61 





587731511532454208 


19.75317 


-0.964191 


20.711 


0.1250 


0.503 


0.120 


0.25 


0.1250 




4963 


7 


4.2 


14.8 


36.5 


14.1 


-14.577 


22.908 


27.428 


1. 


.735 




1.72 





588015507667419343 


19.83900 


-1.181540 


19.716 


1.4350 


0.927 


1.140 


1.53 


1.4350 




4963 


5 


11.1 


43.3 


20.3 


14.1 


-13.737 


26.241 


30.318 


1, 


.565 




1.33 





587727884161581293 


29.96522 


-8.803312 


20.815 


0.1250 


0.597 


0.120 


0.26 


0.1250 




6106 


3 


3.2 


12.0 


32.7 


15.9 


-14.453 


23.032 


27.465 


1, 


.702 




1.77 





587727883893211318 


30.05424 


-8.839227 


20.380 


1.8250 


0.479 


1.440 


2.01 


1.8250 




6106 


2 


5.0 


13.9 


29.1 


15.9 


-14.342 


25.893 


30.304 


1, 


.693 




1.59 





587727178999398623 


30.28336 


-9.408542 


20.874 


2.4850 


0.686 


2.060 


2.76 


2.4850 




3772 


7 


3.5 


11.4 


13.2 


24.2 


-14.278 


26.284 


30.470 


1 


.607 




1.79 





587731512611897674 


32.80721 


-0.133716 


20.592 


1.4450 


0.561 


0.960 


1.66 


1.4450 




2081 


7 


3.2 


10.2 


4.0 


2.7 


-13.872 


26.113 


29.957 


1 


.475 




1.68 





587731513691013244 


45.00239 


0.807781 


16.531 


4-3850 


0.990 


4.190 


4.51 


4.3850 




4145 


7 


0.6 


7.7 


4.0 


7-0 


-13.875 


27.278 


32.749 


2 


.100 




0-05 





587731512083349864 


51.75811 


-0.573625 


20.069 


1-1150 


0.702 


0.620 


1.66 


1.1150 




5810 


2 


8.4 


18.3 


7.8 


6-7 


-13.624 


26.083 


20-515 


1 


.317 




1-83 





587728906008377011 


115.26597 


31.160109 


18.733 


3-8450 


0.971 


3.370 


4.35 


3.8450 




0377 


7 


3.4 


11.7 


26.1 


4-2 


-14.530 


26.488 


31.669 


1 


.989 




0-93 





587728906008376844 


115.31801 


31.184168 


20.169 


4.8350 


0.983 


4.510 


5.52 


4.8350 




0377 


7 


1.8 


9.8 


25.5 


4.2 


-14.595 


26-658 


31-268 


1 


.770 




1.51 





588007005767532860 


116.33380 


39.501996 


20.363 


0.8250 


0.582 


0.390 


1.01 


0.8250 




6111 


1 


5.6 


14.5 


28.7 


8.3 


-14.298 


25.086 


20.585 


1 


.727 




1.59 





588007005767532880 


116.44541 


39.466243 


20.179 


2.0250 


0.738 


1.820 


2.19 


2.0250 




6111 


1 


10.7 


34.7 


35.6 


8.3 


-14.014 


26.333 


30.480 


1 


.592 




1.51 





587731680110052241 


116.73920 


27.665360 


20.778 


4.6050 


0.984 


4.440 


5.19 


4.6050 




3561 


7 


3.2 


10.4 


4.4 


4.7 


-13.833 


27.371 


31.001 


1 


.394 




1.75 





587725470127095818 


118.81215 


41.058484 


19.824 


1.6150 


0.725 


1.440 


1.95 


1.6150 




3032 


2 


11.5 


25.7 


4.1 


8.5 


-13.144 


26.960 


30.331 


1 


.294 




1.37 






(24) 



Note. — The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. (1) SDSS Object ID (2) SDSS R.A. (J2000) (3) SDSS Dec (J2000) (4) SDSS aainh mag.psf i, 
dereddened. (5) Photometric redshift (see Weinstein et al. 2004). (6) Photometric redshift range probability. (7) Lower limit of photometric redshift range. (8) Upper limit of photometric redshift range. (9) Best redshift: spectroscopic 
if different than Zphot- (10) Reference for spectroscopic redshift: S: SDSS, O: ChaMP, R: published reference from NED. (11) Chandra observation ID (obsid) (12) ACIS CCD id (13) Chandra Off-axis angle in arcmin. (14) 99% Counts 
upper limit the 0.3-8keV range (15) Vignetting-corrected exposure time in ksec. (16) Galactic column in units 10^'^cm— 2 (17) log upper limit to the X-ray flux (0.5-8 keV) in erg cm~^ s""'^ (18) log upper limit to the X-ray luminosity 
at 2keV in erg s~^ Hz~^ (19) log optical/UV luminosity t 2500A in erg s~^ Hz~^ (20) aox, the optical/UV to X-ray spectral index. (21) 20cm radio fiux in mjy from FIRST. (22) Radio loudness. (23) Class: 1= BAL, 2= NAL, 3= 
NLSl. (24) Comments 



TABLE 4 
Quasar Sample Univariate Results 



Sample^ 


N 


Mean 


error'' 


Median 


Pmax (%) 




r Distributions'^ 






MainDct 


1135 


1.94 


0.02 


1.93 




RQ 


704 


1.91 


0.02 


1.86 


0.2 


RL 


43 


1.73 


0.05 


1.65 


0.2 


BAL 


15 


1.35 


0.15 


1.30 


0.05 


NAL 


9 


1.68 


0.13 


1.75 


19 


NLSl 


19 


2.01 


0.15 


1.95 


35 


z > 3 


56 


1.80 


0.07 


1.76 


21 






Distributions 






D2L 


1269 


1.421 


0.005 


1.365 




RQ 


680 


1.527 


0.008 


1.452 


0.0 


RL 


31 


1.382 


0.030 


1.392 


0.0 


BAL 


23 


1.717 


0.028 


1.664 


0.0 


NAL 


8 


1.463 


0.056 


1.500 


49 


NLSl 


19 


1.540 


0.080 


1.433 


24 


2 > 3 


47 


1.817 


0.077 


1.786 


55 



^ For each parameter tested, the numeric sample at 
the top is the parent for comparison subsamples be- 
low. The RL subsamples are tested against non-RL 
QSOs from the parent sample. The remaining 3 QSO 
subsamples for each parameter are tested against the 
parent sample excluding all 4 QSO subtypes. 

Error in the mean from the Kaplan-Meier estimator 
as implemented in ASURV. An estimate of the disper- 
sion can be obtained by multiplying this by \/N — 1 

The maximum probability for the null hypothesis 
(of i ndistin guishable samples) from 3 tests described 
in § 16.1.11 Only for Pmax < 5 do we consider the 
distributions significantly different. RL and RQ sam- 
ples are contrasted to each other. Other samples are 
compared to their parent sample (MainDet or D2L)- 
X, where A:=BALs-|-NALs-|-NLS1s, except for BALs, 
whose parent sample is RQ QSOs only. 

These are distributions of "best-PL" measurements, 
best-fit r, which always includes N^"', and also in- 
cludes Nj?"' for 0.5-8 keV counts >200. 



TABLE 5 

Quasar Sample Bivariate Regression Results: 
r{ hkev) OLS 



Sample 


Slope 


error 


Intercept 


error 


MainDet 


-0.1465 


0.0199 


5.8575 


0.5276 


noTDet 


-0.1424 


0.0188 


5.7579 


0.4959 


HiCtNoTRBt 


-0.1213 


0.0497 


5.2473 


1.3333 


NoRBDct 


-0.1528 


0.0262 


6.0343 


0.6980 


hiLoLx 


-0.2336 


0.0388 


8.1937 


1.0337 


HiCt 


-0.1454 


0.0472 


5.8540 


1.2679 



Note. — Samples tested are arranged in the 
same order as in Table [T] OLS refers to the ordi- 
nary least squares regression. 



TABLE 6 

Quasar Sample Bivariate Regression Results: 
hkevi '2500 Bisector 



Sample Slope error Intercept error 



Primary Samples 



Main 


1, 


1171 


0.0170 


-7, 


.5929 


0, 


.6365 


MainDet 


0, 


,9372 


0.0266 


-1, 


.9178 


0, 


.8797 


NoTDet 


0, 


,9719 


0.0259 


-2, 


.9508 


0, 


.8418 


D2L 


1, 


1350 


0.0209 


-8, 


.1162 


0, 


.6329 


D2LN0RB 


1, 


,1667 


0.0238 


-9, 


.0641 


0, 


.7199 


hiLo 


1, 


,2976 


0.0340 


-13, 


.1316 


1, 


.0356 






Other Samples 










NoRB 


1, 


,1502 


0.0146 


-8, 


.5922 


0, 


.6125 


NoRBDet 


0, 


,9359 


0.0253 


-1, 


.8745 


0, 


.8903 


zLxBox 


0, 


,8421 


0.0173 


1, 


.0177 


0, 


.6695 


LoBox 


0, 


,9907 


0.0310 


-3, 


.5191 


0, 


.9357 


zBox 


1, 


,6658 


0.1125 


-24, 


.3936 


3, 


,4355 


zBoxDet 


1, 


,4088 


0.0803 


-16, 


.3512 


2, 


.4564 


D2LSyl 


1, 


,2299 


0.0504 


-10, 


.6975 


1, 


.4746 


D2LN0TRB 


1, 


1908 


0.0196 


-9, 


.7838 


0, 


,5921 


S06 


0, 


,72 


0.01 


4.53 


0, 


.69 



Note. — Samples tested are arranged in the 
same order as in Table [T] OLS refers to the or- 
dinary least squares regression. S06 are results 
from .Steffen et al.. (,2006, ') for comparison. 



TABLE 7 

Quasar Sample Bivariate Regression Results: 



Sample Slope error Intercept error 



Primary Samples 



Main 


0, 


,0598 


0.0066 


-0.2776 


0, 


.1988 


MainDet 


0, 


,0826 


0.0066 


-1.0331 


0, 


.1978 


NoTDet 


0, 


,0732 


0.0089 


-0.7496 


0, 


.2681 


D2L 


0, 


,0610 


0.0085 


-0.3189 


0, 


.2580 


D2LNoRB 


0, 


,0516 


0.0078 


-0.0358 


0, 


.2377 


hiLo 


0, 


1284 


0.0070 


-2.3754 


0, 


.2151 






Otiier Samples 








NoRB 


0, 


,0513 


0.0073 


-0.0239 


0, 


.2217 


NoRBDet 


0, 


,0804 


0.0085 


-0.9667 


0, 


.2546 


D2LNoTRB 


0, 


,0482 


0.0104 


0.0652 


0, 


,3146 


zLxBox 


0, 


,1943 


0.0067 


-4.4484 


0, 


,2033 


LoBox 


0, 


1895 


0.0105 


-4.2956 


0, 


.3176 


zBox 


0, 


1019 


0.0242 


-1.5709 


0, 


.7404 


zBoxDet 


0, 


1119 


0.0267 


-1.9335 


0, 


.8172 


D2LSyl 


0, 


,0058 


0.0175 


1.3253 


0, 


.5160 


S06 


0, 


137 


0.008 


-2.638 





.240 



Note. — Samples tested are arranged in the 
same order as in Table [T] OLS refers to the or- 
dina ry least squar e s regr ession. S06 are results 
from (Steffen et al] (|2006l ) for comparison. 



TABLE 8 

Quasar Sample Bivariate Regression 
Results: r{aox) OLS 



Sample 


Slope 


error 


Intercept 


error 


MainDet 


0.188 


0.106 


1.676 


0.153 


NoTDet 


0.274 


0.108 


1.566 


0.154 


NoRBDet 


0.340 


0.109 


1.466 


0.154 


HiCt 


0.342 


0.202 


1.507 


0.259 


HiCtNoTRB 


0.358 


0.179 


1.529 


0.225 



Note. — Samples tested are arranged in 
the same order as in Table [T] OLS refers to 
the ordinary least squares regression. 



